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ABSTRACT 



Maneuvers employing atmospheric forces to assist in orbital changes hold 
potential for significant fuel savings over purely exoatmospheric propulsive 
methods. The term synergetic has been coined to describe the combination of 
propulsive and atmospheric forces used by a maneuvering flight vehicle. This 
thesis concentrates on non-coplanar synergetic maneuvers using two different 
control methods for various lifting bodies over a range of heating rates and 
orbital speeds. The objective of this thesis is to study the aerocruise and 
aerobang maneuvers. The aerocruise maneuver has been studied for more than 
twenty years and is commonly thought to be the fuel-optimal solution to a 
maneuver flown at a constant heating rate. A new maneuver, the aerobang, has 
recently laid doubt as to the optimality of the aerocruise maneuver. The 
aerobang maneuver demonstrates the ability to yield a higher inclination change 
for a given amount of fuel as compared to the aerocruise maneuver. Within this 
thesis a computer code is developed to model both the aerobang and aerocruise 
maneuvers. It is shown that there exist flight regimes where the aerobang 
method is superior to the aerocruise method. 
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I. INTRODUCTION 



The concept of maneuvering within the atmosphere by flight vehicles moving 
at orbital or near orbital speeds dates back to Sanger's idea of a suborbital 
spacecraft using aerodynamic lift to extend its range [Ref. 1]. While the idea is 
not new, the practical application, using lifting bodies at hypervelocities, may be 
within reach. New materials, propulsion, and engineering techniques have 
revived the idea of a synergetic maneuver. The term synergetic has been coined 
to describe the maneuvering of a spacecraft within the atmosphere using the 
combination of propulsive and aerodynamic forces. A maneuver using 
atmospheric forces to assist in orbital changes holds the potential for significant 
fuel savings over a purely exoatmospheric propulsive method [Ref. 2]. Mission 
applications of these types of maneuvers have been identified in three contexts; 
1) synergetic plane change, 2) planetary missions, and 3) high earth orbit to low 
earth orbit rendezvous [Ref. 3: p. 3]. These three categories encompass the 
totality of useful application of this type of maneuver. This thesis is concerned 
with category one, the synergetic plane change. This application deals with non- 
coplanar orbital transfers, the interest of which is limited to the inclination 
change. A typical flight profile of a spacecraft executing a maneuver of this 
nature would require a deorbit bum which would put the spacecraft into a 
trajectory grazing the atmosphere, an aerodynamic turn, and finally a 
circularization maneuver at the desired altitude (Figure 1). 
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Figure 1. Orbital Transfer Maneuver 
A. THESIS SCOPE 

In the search for the optimum synergetic maneuver, three methods have been 
studied: the aeroglide, the aerocruise, and the aerobang. The aeroglide 

maneuver uses no propulsive force during the atmospheric interaction, 
theoretically giving this method an infinite inclination change per kilogram of 
fuel used. The penalty for the aeroglide maneuver is the extremely high heating 
rates caused by the dense atmosphere from the required low altitudes. The 
aerocruise maneuver was developed to overcome the heating rate problem 
created by the aeroglide method. In order to maintain a constant heating rate the 
aerocruise maneuver uses a combination of propulsion and aerodynamics to 
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achieve an inclination change. For more than twenty years, this method has been 
studied to determine the optimal regimes of operation for this maneuver. A 
recent study has raised doubt to the optimality of the aerocruise maneuver and 
proposed another maneuver based on a maximum throttle setting and a 
modulation of the angle of attack [Ref. 3]. This new aerobang maneuver 
demonstrates the ability to yield a higher inclination change per kilogram of fuel 
used when compared to the aerocruise maneuver under like conditions. 

This thesis analyzes two control methods of synergetic maneuvers, the 
aerocruise and the aerobang method. Each of the two maneuvers have similar 
goals of maintaining a constant heating rate while executing the aerodynamic 
turn. An evaluation, analysis, and comparison of the two maneuvers over a 
range of maximum lift to drag ratios ((L/D)max) and heating rates will shed light 
on the capabilities and inadequacies of each control method. The analysis is 
limited to the actual maneuver itself at the necessary altitude required for its 
execution. It was assumed that the the method of entry or exit from the 
atmosphere for the two maneuvers would be identical and would not be a factor 
in the overall efficiency of the maneuver. 

The study includes a program developed specifically to model the aerocruise 
and aerobang control laws. The performance of each of these strategies are 
analyzed over a variety of orbital speeds and altitudes. The intent is to develop 
insight into the range of capabilities for each method and conditions under which 
aerobang might be more efficient than aerocruise. 

B. BACKGROUND 

Synergetic non-coplanar orbital transfers utilize a portion or all of the lift 
generated by the body in a direction perpendicular to the orbital plane. A force 
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applied in this direction will cause a change in the inclination of the orbit. The 
effect of the force on the inclination is maximized at the ascending and 
descending nodes of the orbits but has no effect at the apex of the orbit (i.e, 90 
degrees from the nodal crossing). This maximizing effect at the node 
necessitates that the maneuver be completed within a short time span relative to 
the orbital period of the flight vehicle. Three methods of inclination change 
have been proposed so far: aeroglide, aerocruise, and aerobang. While this 
thesis deals with only two of these, the third method, aeroglide, would be of use 
in bringing the spacecraft to the appropriate position in which to commence the 
maneuver. 

1. Aeroglide 

During a non-coplanar transfer the aeroglide maneuver seeks to change 
the orbital inclination of the spacecraft. The maneuver consists of three parts, 
the deorbit, the aerodynamic turn, and the reorbit. The first and last segments of 
the maneuver are purely propulsive and generic to any of the three methods 
discussed within this thesis. To describe the aerodynamic turn of the aeroglide 
maneuver, it must be understood that no propulsive forces are used within the 
atmosphere. In order to maximize the plane change, the maneuver is most 
efficient at an angle of attack that produces the maximum lift-to-drag ratio. The 
lifting body is oriented with the generated lift vector normal to the spacecraft's 
orbit. To generate the necessary lift for an effective plane change, the vehicle is 
flown at an orbital speed and a low altitude where the atmosphere is very dense. 
This flight profile results in an extremely high heating rate for the flight vehicle, 
possibly beyond the capabilities of the spacecraft. If the heating rate problem 
can be overcome, the aeroglide maneuver requires fuel expenditures only to 
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deorbit and reorbit, significantly reducing the fuel required to effect the same 
plane-change using a pure propulsive method. [Ref. 5: p 104] 

The extreme heating rate required by the aeroglide maneuver hinders its 
use as an effective method for non-coplanar transfers. The heating rate profile 
of the maneuver begins at a low value as it enters the atmosphere and builds to its 
highest value at the minimum orbital altitude, A possible use for the aeroglide 
maneuver would be in conjunction with another synergetic plane change 
maneuver. A proposed flight profile would be: after the deorbit bum the 
spacecraft would assume an aeroglide profile until the predetermined altitude, 
velocity, and heating rate were achieved, then the follow-on maneuver would be 
executed, and finally an aeroglide maneuver to exit the atmosphere and reboost. 
Both aeroglide portions of the flight path would not be at the optimum glide 
profiles for an aeroglide maneuver, but some gain of inclination change would 
be realized until the selected conditions for the aerocruise or aerobang maneuver 
are reached, 

2. Aerocruise 

The aerocruise synergetic maneuver couples the aeroglide maneuver 
with a propulsive force during the atmospheric interaction. This maneuver is 
made up of three parts, the deorbit, the aerodynamic turn, and the reorbit, 
similar to the aeroglide maneuver, the difference being in the method used 
during the aerodynamic turn. A typical flight profile requires the vehicle to 
remain on the heating rate boundary; this dictates the aerocruise maneuver to be 
flown at a specific altitude and speed. Graphically, this profile appears as a point 
on a speed versus radius graph (Figure 2). This point would just touch the 
maximum allowable heating rate for the spacecraft. The altitude and speed are 
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maintained by a combination of engine throttling and modulation of the flight 
vehicle angle of bank. The speed is held constant by setting the thrust equal to 
the drag encountered at the flight altitude. Altitude is maintained by banking the 
spacecraft to allow for a component of the lift force generated to counteract the 
apparent centrifugal force (Figure 3). This method allows for a plane-change 
maneuver without incurring the heating rate problems of the aeroglide 
maneuver. Conventional wisdom requires that the angle of attack be set at 
(L/D)max in order to maximize the efficiency of the maneuver [Ref. 3: p. 5]; 
however, recent studies [Ref. 6: p. 517] have shown that the most efficient (i.e. 
largest inclination change for a given expenditure of fuel) maneuver is not at 
(L/D)max, but at a lower L/D and a higher angle of attack. 



6 




Figure 3. Graphical Representation of Aerocruise Maneuver 
3. Aerobang 

This method of control is similar to the aerocruise model in that the 



heating rate of the vehicle is to remain constant, but this is where the similarity 
ends. The aerobang flight profile allows the vehicle to change its altitude and 
speed as necessitated by the control law governing the maneuver (Figure 2). As 
mentioned previously, the most effective use of the lifting force for the orbital 
plane change is to orient the force normal to the orbital plane; this dictates an 
angle of bank of 90 degrees for the spacecraft. If the angle of bank is to remain 
fixed, the only variables left for modulation are the thrust of the vehicle and the 
angle of attack. The aerobang maneuver modulates the angle of attack and fixes 
the thrust at the maximum available value. It can be shown that for certain 
conditions, the aerobang maneuver will achieve a higher inclination change for a 
given fuel expenditure than the aerocruise maneuver. 
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II. FLIGHT VEHICLES 



To properly understand the similarities and differences between the 
aerocruise and aerobang maneuvers, a parameter analysis of more than one flight 
vehicle would be necessary. Ideally, the flight vehicles would be similar in every 
respect except one, and changes observed on the data output could then be 
attributed to that differing parameter. Early investigators of synergetic 
maneuvers came to two main conclusions: first, the lifting body should operate 
at (L/D)max. and second, the L/D should be greater than one to be superior to an 
all propulsive maneuver [Ref. 3: p. 5]. The emphasis on (L/D)max brought this 
investigation to select two vehicles for comparison, one which would operate at 
(L/D)max of approximately 1,8 and the other at 2.6. Flight vehicles which have 
been designed and have data available provide a sense of reality to the analysis, 
where as an imaginary flight vehicle and data might skew the study. The two 
flight vehicles investigated in this thesis are the ERV (Entry Research Vehicle) 
and the MRRV (Maneuverable Reentry Research Vehicle). 

A. ENTRY RESEARCH VEHICLE 

The ERV was designed with the idea of investigating and exploiting 
maneuvers involving long downrange, high crossrange, and synergetic plane 
changes. The vehicle is 7.62 m long and has a wing span of only 3.96 m, 
allowing the vehicle to be launched from the cargo bay of the space shuttle. The 
slenderness ratio of the fuselage is approximately 0.167. Figure 4 shows a three- 
view sketch of the ERV. Three Marquardt R-40-B rocket motors make up the 
propulsion system, providing a total of 14679 Newtons of thrust at a specific 
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Figure 4. Three-view Diagram of ERV 
impulse (ISP) of 295 seconds. Separation weight from the Space Shuttle is 
planed for 5443 kg with 50% of that made up of fuel. [Ref. 7: pp. 489-90] 

For the purpose of modeling within a computer program, the aerodynamic 
data were approximated by a least squares binomial fit (Figures 5 and 6). Wind- 
tunnel data for the ERV were obtained from the Martin Marietta Corporation 
[Ref. 8] and included coefficient of lift and drag data up to Mach 10. These data 
were used for flight velocities from 8.3 km/sec down to 7.1 km/sec. Although 
these velocities correspond to Mach numbers of 26 to 22 at the maneuvering 
altitude, the data from Mach 10 can realistically apply because of the limiting 
characteristics exhibited by the bow shock and the pressure distribution over the 
lifting body [Ref. 9: pp. 387-388]. The curve fit was done using a linear 
statistical model [Ref. 10] and then the coefficients were modified to obtain an 
L/D vs angle of attack similar to the actual data taken from the wind tunnel tests. 
From the L/D curve (Figure 7), we see that (L/D)max occurs at approx 10.5 
degrees with a value of 1 .8. 
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Figure 5. Coefficient of Lift vs Angle of Attack for ERV 




Figure 6. Coelficient of Drag vs Angle of Attack for ERV 
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Figure 7. Cl/Cd vs Angle of Attack for ERV 



B. MANEUVERABLE REENTRY RESEARCH VEHICLE 

Although the MRRV is similar in size and weight to the ERV, that is where 
the similarity ends. The MRRV has a maximum lift to drag ratio on the order of 
2.6, much higher than the ERV. Figure 8 show a three-view sketch of the 
MRRV and when compared to Figure 4 the differences between the MRRV and 
the ERV become apparent. The MRRV is a much more slender design; the 
slenderness ratio of the fuselage is approximately one eighth. The overall length 
is 7.62 m and because of the sharper design of the vehicle it has an effective wing 
area of 11.7 square meters, 4.7 square meters less than the ERV. The propulsion 
system is identical to the ERV using three Marquardt R-40-B rocket motors, 
providing a total of 14679 Newtons of thrust at a specific impulse of 295 
seconds. [Ref. 3: p. 7] 
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Figure 8. Three-view Diagram of MRRV 
The method used to find the coefficient of lift and drag curves for the ERV 
was employed to find those for the MRRV. Figure 9 shows a plot of the curve 
fit for the coefficient of lift with the wind tunnel data obtained for the MRRV, 
and Figure 10 displays the coefficient of drag. It is interesting to note that even 
though Figure 1 1 shows the MRRV to have a higher lift to drag ratio than the 
ERV, a comparison of the actual coefficient of lift and drag show the ERV to 
have higher actual values. The lower lift and drag coefficients coupled with the 
decreased wing area of the MRRV will be of importance during the analysis to 
follow. 
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Curve fit ^ Data 



Figure 9. Coefficient of Lift vs Angle of Attack for MRRV 




Curve fit O Data 



Figure 10. Coefficient of Drag vs Angle of Attack for MRRV 



13 




Figure 11. Cl/Cd vs Angle of Attack for MRRV 
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III. EQUATIONS OF MOTION 



The equations of motion for a hypervelocity lifting body are derived, free of 
singularities. Six specific variables are utilized to describe the spacecraft's 
position and velocity. Referenced to an inertial frame XYZ, the vehicle's 
position is represented by spherical coordinates r (radius), 0 (right ascension), 
and <j) (declination) and the vehicle's motion as speed V, and direction by the 
angles y (flight path angle) and \j/ (heading angle). The flight path angle, y, is 
positive above the local horizon and negative below. The heading angle is the 
angle the horizontal projection of the velocity vector makes with the local 
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latitude. Figure 12 graphically shows the six variables and their relationship to 
each other. 

Assumptions made in the derivation of these equations are: a spherical, non- 
rotating, Earth and negligible atmospheric winds. The Earth's oblateness is not 
considered due to the short nature of the atmospheric interaction by the 
spacecraft. Future studies may find a need to include these effects. Additionally, 
the thrust vector of the spacecraft was assumed to be in line with its longitudinal 
axis. 

A. DERIVATION OF EQUATIONS 

Three coordinate systems are defined for the derivation of the equations of 
motion. The inertially fixed system, XYZ. A rotating system, xyz, in which the 
origin is coincident with that of the inertial system and the positive x axis 
remains fixed to the spacecraft. Finally, a system x'y'z' which is fixed to the 
spacecraft with the y’ axis oriented along the velocity vector and the x' axis 
within the xy plane, the origin of the system lies at the end of the radius vector in 
the xyz system, (Figure 12). The transformation from the inertial system into 
the xyz system is done using a body 3-2 rotation through angles 0 and (() 
respectively. Additionally, to obtain the xyz system from the x'y'z’ a body 3-1 
rotation must be done using angles y and y. Using the following DCMs 
(direction cosine matrix) all of the parameters can be described in the rotating 
xyz system. 

The DCM for the transformation from the inertial system to xyz system is 
expressed as the multiplication of two matrices the right matrix representing the 
body-3 rotation and the left matrix representing the body-2 rotation 
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-sin (]) 


0 


cos (]) 


- 0 


0 


1 - 


-Z- 



Consolidating the two matrix transformations for a body 3-2 rotation, we have 



*x 




cos 0 cos ([) 


sin 0 cos ([) 


sin ({) 


'X' 


y 


- 


-sin 0 


cos 0 


0 


Y 


LzJ 




-cos 0 sin (j) 


-sin 0 sin <]) 


cos ({) 


-Z- 



Similarly, the DCM for the transformation from the x’y'z' spacecraft system to 
xyz system is expressed below: 



(3.3) 



X 




1 


0 


0 ■ 


’ cos Y sin Y 


0 ■ 


x'‘ 


y 


= 


0 


cos \|T 


-sin Y 


-sin Ycos Y 


0 


/ 

y 


.zi 




0 


sin Y 


COSY 


0 0 


1 


z' 



Completing the algebra the body 3-1 rotation, we get. 



X 

y 

LZI 



COS Y sin Y 0 

-sin Y COS \|T cos y cos \|r -sin \\f 

-sin Y sin \|r cos y sin \|r cos \\f 



X 

/ 

y 

z 



(3.4) 



For a system of equations to completely describe the spacecraft's motion and 
position at any moment in time, an equation of the time derivative for each of the 



six variables must be established. By equating the spacecraft’s velocity vector to 
the time derivative of the spacecraft’s position vector, the derivatives for r, 0, 



and (j) can be expressed as functions of r, 0, (|), V, y, and \}/. The time derivatives 
for V, Y, and \|/ are also found by taking the time derivative of the spacecraft’s 



velocity vector and equating them to the forces acting on the spacecraft using the 
relationship 
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The vectors for the spacecraft’s position and velocity in the xyz system are 



r = r 1 



(3.5) 



and 

A A 

Vj' = V sin Y i + V cos y cos V|/ j + V cos y sin \|/ k 
Additionally, the angular velocity of the spacecraft is 

= sin 4> ^ i - ^ j + cos (}) ^ k 
^ dt dt ■’ ^ dt 

Taking the derivative of Equation 3.5, yields an equivalent velocity vector, 

V = ^ =^i + ( Qxf ) = ^i + r cos (}) + r^k 

dtxYZ dt dt dt dt 



(3.6) 



(3.7) 



(3.8) 



Equating the components of Equations 3.6 and 3.8 gives the first three equations 
of motion. These three equations describe the spacecraft's position in spherical 
coordinates for any time 

(3.9) 



^ = V sin y 
dt 



V cosy cosy 



dt 



r cos <() 



d(J) _ V cos y sin y 
dT~ r 



(3.10) 



(3.11) 



Similarly taking the derivative of the velocity vector (Equation 3.6) the final 

three equations of motion can be found. The last three equations will describe 

the spacecraft's movement through space at any particular time. 

d(Vsiny)* d( Vcosycosvi/) * 

' dt 



dt 



dt 



dt 



dj d ( V cos y sin \ii) ,, . au 

+ V cosy cos H ^k + V cosy sin \i/^ 

^ dt dt ' ^ dt 
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Expanding, we have 



dV = 
dt 



sin y^ + W cos 
' dt ' dt. 



1 + 



V sin Y ( cos <}) ^ j + ^ k ) 
' ^ dt •' dt . 



jv dy . d\if 

cos y cos \|/ ^ - V cos \it sin y — - V cos y sin \ii — — j 
' ^ dt ^ dt ' ^ dt . •* 



V cos Y cos XI/ ( -cos (b ^ i -i- sin 6 ^ k) 
' ^ ^ dt dt 



dV • • dy d\|/ 

cos y sin \i/ - V sin w sin y — + V cos y cos \i/ 

' ^ dt ^ ' dt ' ^ dt J 



d(J) A HA * 

V cos Y sin \i/ ( - — 1 - sin d) 1) 
' ^ dt dt 



(3.12) 



Replacing d9/dt and d(})/dt in Equation 3.12 using Equations 3.10 and 3.1 1 yields 



dV^ 

dt 



HV dy cos^ y cos^ w cos^ y sin^ \ii 

sin Y^ + V cos Y^ ‘ — 

' dt ' dt r r 



Hv dy . dw 

+ [ cos Y cos \i/ ^ - V cos \i/ sin Y — - V cos y sin \it -r- 
' ^ dt ^ ' dt ' ^ dt 

sin y cos y cos \|t V cos^ y sin y cos y tan ()) ^ ^ 

r r -I J 



-t- [ cos Y sin y - V sin y sin y ^ -i- V cos y cos y ^ 

' ^ dt ^ ' dt ' ^ dt 

sin y cos y sin y cos^ y cos^ y tan d) , 

-I- — — ^-1- — — -] k 



dy 



(3.13) 

r r 

The forces acting on the spacecraft in the x'y’z' coordinate system can be 

expressed as (Figure 13) 



F^'=( T sin a + L ) cos a 

Fy'= T cos a - D 
Fj'= ( T sin a -I- L ) sin a 



(3.14) 

(3.15) 

(3.16) 
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Restating the fundamental equation of motion as a derivative of velocity with 
respect to time 



F = niai = ma^ 



(3.17) 




Figure 13. Spacecraft Forces During Atmospheric Interaction 
Using the DCM from Equation 3.3 the forces from the spacecraft's frame of 
reference (x'y’z’) can be transformed into the rotating frame of reference (xyz) 
using 



Fx 

Fy 

Fz 



cos Y sin Y 

- sin Y cos \|/ cos Y cos \|/ 

- sin Y sin \(/ cos Y sin \\f 



0 

- sin Y 
cos 




(3.18) 



With the above transformation matrix and using Equations 3.14, 3.15, and 3.16, 



the forces acting on the spacecraft can be expressed in the rotating system as 



= ( T sin a + L ) cos a cos Y + ( T cos a - D ) sin y * nia g 



(3.19) 



Fy = - ( T sin a + L ) cos a sin ycos \(/ 



+ ( T cos a - D ) cos ycos y - ( T sin a + L ) sin a sin y 



(3.20) 



20 



and 



= - ( T sin a + L ) cos o sin y sin \/ 

+ ( T cos a - D ) cos Y sin \j/ - ( T sin a + L ) sin a cos y (3 21) 

Now setting the components of Equation 3.17 equal to the spacecraft force 
equations expressed in the rotating coordinate system, three independent 
equations are obtained for evaluation and reduction namely 

A 

i : ( T sin a + L ) COS a COS Y + ( T COS a - D ) sin Y - m g = 



m„ 



dV ^7 ^7 V^cos^Y 

sin Y ^ + V cos Y 

dt dt r 



( 3 . 22 ) 



j : - ( T sin a + L ) cos a sin Ycos y + ( T cos a - D ) cos Ycos 

- ( T sin a + L ) sin G sin y = [ cos ycos y ^ - V cos y sin y^ 

dy \/2 

- V cos Y sin y cos y cos y( sin y - cos y sin y tan (J) ) ] 

k:-(Tsina + L) cos a sin y sin y + ( T cos a - D ) cos y sin y 

+ ( T sin a + L ) sin G cos y = [ cos y sin y ^ - V sin y sin y^ 

+ V cos Y cos y ^ ^ cos y ( sin y sin y+ cos y cos^ y tan <]) ) ] 



(3.23) 



(3.24) 



In order to obtain the last three equations of motion. Equations 3.22, 3.23, and 
3.24 must be algebraically manipulated and reduced. Taking Equation 3.23 and 
multiplying by cos(y) and Equation 3.24 and multiplying by sin(y), then adding 
the two resulting equations together gives 



( T cos a - D ) - g sin Y = ^ 



(3.25) 
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Substituting this result back into Equation 3.22 and reducing yields 

, , dy ( T sin a + L) cos a - m g cos y (3 26 ) 

V — = h COS 7 ^ 

dt r ' 

A final substitution into Equation 3.23 gives us the sixth equation of motion for 

the flight vehicle, namely 



dy ( T sin a + L) sin a . 

V — ^ cos y cos w tan (b 

dt m cos y ^ 



(3.27) 



Making one final reduction by defining variables Ar, As, and Aw which are 
the components of the acceleration of the flight vehicle in the radial, tangential, 
and binormal directions respectively, we have from Equations 3.14, 3.15, and 
3.16, 



and 



^ _ ( T sin a + L) cos a 

^ _ T cos g - D 
S m 



Aw = 



( T sin a + L) sin a 






The last three equations of motion can now be written as 



dV 

dt 



= As - g sin y 



V^ = AR + (^-g)cosy 
and 



ydy_ Aw 
dt cosy 



+ 






cos y cos tan (]) 



(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 
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B. DERIVATION OF THE INCLINATION 

The six equations of motion. Equations 3.9-3.11 and 3.31-3.33, when 
integrated numerically yield the spacecraft’s position in spherical coordinates at 
any time during the orbit. For facilitating a comparison between orbits, a 
relationship between the orbital elements and the spherical coordinates is 
necessary. In non-coplanar transfers, an orbital element which indicates 
performance of the maneuver is inclination. The necessity to convert from 
spherical coordinates to inclination of the orbit for the purpose of performance 
evaluation is the basis for this derivation. As shown in Equations 3.5 and 3.6 the 
radius and velocity vectors can be written in terms of the rotating frame of 
reference. The angular momentum vector ( h ) is defined as the cross product of 
the radius and the velocity vector; thus 

— -♦ A y«S, 

h = r X V = - V r COS Y sin \j/ j + V r COS Y COS \|r k (3.34) 

The magnitude of the momentum vector is, consequently 

|h] = V( V r cos^ Y sin^ \|/ + ( V r cos^ Y cos^ \|/ (3.35) 

giving 

|h| = VrcosY (3.36) 

The inclination of an orbit is defined as the the angle between the angular 
momentum vector and k, the direction of the Z axis of the inertial coordinate 
system. The inclination can be found by defining k within the rotating xyz 
system and then using the dot product definition to find the angle between the h 
and the k vectors. Since 

— * A y-s. 

k = sin <[) i + cos ({) k (3.37) 
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h-k = 



cos 1 



(3.38) 



and 

V r cos Y cos y cos <|) cos i (3.39) 

Utilizing Equations 3.36 and 3.39, the cosine of the inclination can be found 
explicitly in terms of y and (j), thus 

cos i = cos y cos (|) (3-40) 

This result, exhibits no constraint on whether the path of the orbiting vehicle is 
Keplerian or not. 
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IV. CONTROL LAWS 



Aerocruise and aerobang each have unique methods for controlling the space 
vehicle's trajectory during an atmospheric encounter. Aerocruise uses two 
distinct rules to maintain the proper flight path. First, the spacecraft’s velocity is 
controlled by setting the component of thrust in the direction of travel to exactly 
counteract the vehicle's drag. Second, the altitude is maintained constant by 
modulating the flight vehicle's angle of bank. This is done by balancing the 
vertical component of the lift and thrust to the vehicle's weight minus the 
centrifugal force. The method of control used by the aerobang maneuver is less 
complicated, in that only one parameter of the flight vehicle is modulated during 
the maneuver. The angle of bank during the aerobang maneuver is set at a 
constant value, typically 90 degrees; this maximizes the inclination change caused 
by the aerodynamic lift. Additionally, the thrust is not modulated but set to the 
maximum available thrust. The modulated parameter is the flight vehicle's angle 
of attack, and the spacecraft is not constrained by a particular speed or altitude, 
but to the maximum allowable heating rate. 

A. AEROCRUISE CONTROL LAW 

Early investigations of the aerocruise maneuver dictated that superior results 
could be achieved if the spacecraft operated at (L/D)max [Ref. 3: p. 5]. The 
development of the control law for the angle of bank required to maintain a 
constant altitude and velocity, and hence a constant heating rate, utilizes this 
concept. Counterbalancing the drag on the flight vehicle with the thrust 
produced demands a prediction of the drag encountered at various altitudes. The 
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density in the U.S. Standard Atmosphere, 1976, [Ref. 11: p. 7] is approximated 
by an exponential behavior in the altitude range from 50 km to 120 km, namely 

p = Ps (4.1) 

where 



ps = 3.0968x10-^ kg/m3 
p= 1.41x10-4 m-i 
rj = 6438 km 



This density approximation coupled with the speed of the spacecraft and the 
coefficient of drag on the vehicle, at (L/D)max» yields a drag force which is 
constant throughout the aerocruise maneuver. Looking back at Figure 13, the 
thrust required for aerocruise to balance the drag force by the component of 
thrust along the spacecraft’s direction of motion is given by 



where 



T — D 
cos a 



(4.2) 



D = ipV"sCD (4.3) 

In order to calculate the required angle of bank for the aerocruise maneuver it is 
necessary to look at Equations 3.28 and 3.32, which are repeated here for 
convenience 



^ _ ( T sin a + L) cos a 

Ar- is;; 

V^ = AR + (^-g)cosy 



(3.28) 

(3.32) 



The flight path of the vehicle during the aerocruise maneuver is circular and thus 
the flight path angle (y) is zero and will remain zero throughout the maneuver. 
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Combining this information, with the substitution of Equation 3.28 into Equation 
3.32 and solving for the angle of bank (ct) yields 



a = cos 






(4.4) 



( T sin a + L) 

Obviously, from Equation 4.4, the angle of bank becomes 90 degrees when 
the spacecraft's velocity is equal to the velocity required to maintain a circular 
orbit at the given altitude. As the speed of the spacecraft becomes higher than 
circular speed, the angle of bank also increases beyond 90 degrees, causing the 
lift vector to be directed Earthward to maintain the circular path of the vehicle. 
Likewise, as the spacecraft's velocity decreases below circular speed, the angle of 
bank decreases and the lift vector is directed upward to sustain the circular path. 
Examining Equation 4.4, it is clear that there is a limitation in the angle of bank 
which can be used to control the flight path; this limitation occurs when the 
absolute value of the terms within the square bracket exceeds 1. In this case, 
there is not enough lift generated to maintain the desired flight path. To correct 
for this condition, the vehicle would have to abandon flight at (L/D)max 
choose a higher angle of attack in order to generate the necessary lift. 



B. AEROBANG CONTROL LAW 

Differing from the aerocruise method, the aerobang maneuver allows the 
vehicle to vary the altitude and velocity in order to maintain a constant heating 
rate on the spacecraft. In order to determine how the angle of attack is 
modulated to keep the flight vehicle at a constant heating rate, the modeling of 
the heating rate must first be discussed. Heating rate is modeled as a function of 
atmospheric density (p) and vehicle speed (V), by the expression [Ref. 9: p. 291], 
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0 = k p" V>" 



(4.5) 



where K, n, and m are constants. 

Coupling Equation 4.5 with Equation 4.1, the exponential atmospheric 
model, the speed of the spacecraft can be shown to be a function of the vehicle’s 
geocentric radius for any given heating rate. Equation 4.6 shows this 
relationship and is the basis for the aerobang control law. 



V = 



Oo 



-iJL 

m 



(4.6) 



k pSe-Pn(r-r») 

The derivative with respect to time of Equation 4.6 is necessary in order to 
couple the equations of motion for the spacecraft to the constant heating rate 
flight path. Thus, 



dV = 
dt 



' Qo ‘ 


m 


pn 


k pS. 


Cm" 


L m J 



P n t 
-!_r ^ 

® dt 



dV = 

dt 



pn 



m J 



vdl 

dt 



(4.7) 



(4.8) 



Combining Equation 3.9 with Equation 4.8, dr/dt can be eliminated and dV/dt 
can be expressed as a function of V and y 

pn 



dV _ 



sin 7 



(4.9) 



dt >■ m 

Substituting Equation 3.29 into Equation 3.31 and eliminating the variable As, 
another equation for dV/dt is derived (Equation 4.10). The two Equations, 4.9 
and 4.10, yield a final relationship 



dV _ T cos g - D . g sin -Y 
dt m. ^ ' 



(4.10) 



Equations 4.9 and 4.10 together enable the evaluation of the angle of attack 
for any given flight condition, from the equation 
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(4.11) 



Bn 2 

0 = T cos a - D * nig sin Y ( g + V ) s f ( a ) 

However, since the drag term, D, is also a function of a, solving for a is not as 
straightforward as it might at first appear. 

1. Mechanics of Convergence 

Equation 4.11 is made up of three distinct terms, a thrust component, a 
drag component, and a component which depends on the vehicle’s velocity and 
flight path angle. The first two of these terms are dependent on the angle of 
attack of the vehicle, while the third term remains constant for a given flight 
condition, independent of the angle of attack. Figure 14 graphically show the 
relationship of these terms in solving for the angle of attack of the flight vehicle. 
The terms which contain the flight path angle and the drag are summed and 




Figure 14. Graphical Example of Angle of Attack Solution 
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plotted together with the thrust term, the intersection of these two lines occur at 
the proper angle of attack to maintain a flight path along a constant heating rate 
line. 

For a given flight vehicle, the thrust curve shown in Figure 14 is fixed 
for any flight condition, while the drag and flight path terms increase with 
increasing flight speed. As the flight speed continues to increase, the converged 
angle of attack needed to maintain a stable flight path decreases, and at some 
point there exists a flight condition where the two curves would not intersect. At 
this point, the thrust curve lies below the curve representing the sum of the drag 
and flight path terms, and no solution exists for an angle of attack. 

2. Newton’s Method of Convergence 

The nature of Equation 4.11 does not allow a direct solution for the 
angle of attack, and another method must be used. An iterative process 
developed by Sir Isaac Newton lends itself well to applications involving 
computer methods [Ref. 12: p. 142]. The method involves following the tangent 
of the function down to the zero of the function as expressed by 

/'(«„)’=0 (4.12) 

where, 

r(a) = -Tsina-^ = -Tsina-(C 5 + 2C6|o|) jP V^S (4 13) 

and C 5 and C 6 are the coefficients of the first order and second order terms in 
the binomial curve fit of drag vs. angle of attack. The iterative nature of this 
process requires an accurate initial guess for the angle of attack in order to 
reduce the time of computation. 
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3. Initial Guess of Angle of Attack 

Properly setting the initial angle of attack close to the correct value can 
significantly decrease the required number of iterations needed to converge to 
the solution. The method employed here involves simplifying the information 
already presented in order to obtain a direct solution. Two assumptions are 
made in order to provide this starting point for the solution of a. These two 
assumptions are that initially the acceleration component in the tangential 
direction is nearly zero and that the angle of attack is small enough that its cosine 
is approximately unity. Applying these two simplifications in Equation 3,29 
leaves the expression that thrust (T) equals the drag (D), Therefore 

ipV^S ' •’ 

yielding the solution 



C5 + 



ci + 4 



a = 



ipv^s 



-C4 



C6 



2Ce 



(4.15) 



While this is only an approximation, it enables the initial angle of attack 
to start closer to the solution than if some arbitrary value for a were selected. 
As will be shown later, these assumptions are close to the actual value for a flight 
vehicle with a circular or near circular velocity. 
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V. PROGRAM DEVELOPMENT 



A computer program was developed using VAX FORTRAN to simulate a 
lifting body within the atmosphere. The program utilized a fourth order Runge- 
Kutta routine to integrate the differential equations of motion derived in Chapter 
III. Equations of motion for the simulation were left in a dimensional form 
because the accuracy of the computations was considered adequate using double 
precision variables within the program. The limited gain in accuracy and 
computational time using non-dimensional equations and programing techniques 
was considered insufficient to justify the increased time for implementing the 
equations in a non-dimensional form. Redundant checks were made both internal 
and external to the program to ensure the integration time step of the Runge- 
Kutta routine was kept small enough to insure the accuracy of the program. 

The program was written in a modular format with the main program 
controlling six subroutines. Three of the subroutines, CNTRL, ACEL, and 
ORB, deal with the flight mechanics of the spacecraft. Two of the subroutines, 
HDR and WRT, control the output of the program and the last subroutine, RK4, 
is a standard fourth order Runge-Kutta routine. The source code of the program 
is given in Appendix A. 

A. MAIN PROGRAM 

The main program, pages 68 through 70 of Appendix A, acts as a control 
center for the subroutines, and computations within the main program are kept 
to a minimum. Only the calculation of the spacecraft mass is done within 
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Figure 15. Flow Chart for Main Program 



33 



the main program. Figure 15 is a flow chart of the main program operations. 
Two input files were used in order to separate the orbital parameters from the 
vehicle aerodynamic parameters. Sample input files are included in Appendix A. 
One output file was used for both a hard copy text and input into a separate 
graphics program. A quick health check of the program was displayed to the 
terminal for verification of proper inputs and progress through the program. 
The first call to subroutine CNTRL and ACEL are for output purposes only, the 
values at time zero are displayed in the output file after calling these two 
routines. The main loop of the program proceeds through the aerodynamic, 
atmospheric, and orbital dynamic models at each integration time step, the time 
interval is then compared with the time print interval and an output written if 
necessary. The last part of the loop is to change the spacecraft mass in 
accordance with the time interval, thrust on the spacecraft, and the specific 
impulse of the rockets. A double check is used to determine the end of the 
program: if either the final time or final mass is reached, the program is 

terminated. 

B. SUBROUTINES 

1. Subroutine CNTRL 

This subroutine, as the acronym implies, deals with the control laws used 
for the different atmospheric interaction models. A flow chart of this subroutine 
is shown in Figure 16. Both the aerobang and aerocruise models were included 
in the subroutine and each was accessed through the variable CSE. The 
atmospheric model, which is common to both models, is executed prior to the 
control functions. A test is then done in order to determine which model to turn 
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Figure 16. Flow Chart of CNTRL Subroutine 
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on; for convenience, the order of discussion will be the same as their appearance 
within the subroutine. 

a. Aerobang Control Model 

This model is made up three distinct parts. The first part calculates 
an initial guess of the angle of attack for the flight vehicle. While the 
convergence routine will converge on the correct angle of attack using any initial 
guess, a good first guess is determined to decrease the computation time of the 
program. The second portion of the convergence routine uses Newton’s method 
of convergence. First, the coefficients of drag and lift for the flight vehicle for 
a given angle of attack are computed. This information is then used to calculate 
the value of the f(a) in. Equation 4.11, and its derivative using. Equation 4.13. 
These two values dictate the change in the angle of attack and the whole process 
is repeated until there is little change to the angle of attack in sequential loops. 
The final part of the model is an automatic internal check. It is quite possible for 
the model to have no angle of attack on which to converge. Two possibilities 
exist for this happening: one, the model is unable to maintain a constant heating 
rate for the given flight conditions and two, that the integration step size is too 
large, forcing the vehicle off the constant heating rate curve. If there were no 
check within the routine, the output would be questionable, or, the routine might 
enter into an infinite loop looking for a nonexistent answer. 

b. Aerocruise Control Model 

The second model within the CNTRL subroutine incorporates 
Equations 4.3 and 4.4. The model sets the angle of attack for (L/D)max and lift 
and drag are calculated. Using this drag prediction, the thrust is set to a value 
which will exactly match the drag. Using all these values, the angle of bank 
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necessary for circular flight is found. A check is also done to assure that there is 
an angle of bank which will satisfy the given flight conditions; if one does not 
exist, then the program is stopped at this point . 

2. Subroutine ACEL 

This simple routine computes the accelerations on the flight vehicle, 
excluding gravity, using Equations 3.28-30. These equations have been isolated 
within their own subroutine for future program users. Accelerations 
experienced by the flight vehicle during a maneuver can be evaluated using the 
information provided by this routine. 

3. Subroutine ORB 

Subroutine ORB calculates the differential values contained within the 
equations of motion derived in Chapter III. These equations are defined within 
the subroutine using the dimensioned variables X and XDOT. This is done in 
order to make it easier to pass the variables and their derivatives to the Runge- 
Kutta routine and to make the routine more generic in nature. 

4. Subroutines HDR and WRT 

These are the output subroutines and HDR is called only once during 
each execution of the program. The subroutine HDR prints the header at the top 
of the output file. This header contains some of the selected input data in order 
that one may readily associate the conditions that created it. The subroutine 
WRT prints the output data and is called at every time print interval as specified 
by the user. It is necessary to limit printing of the output to every 1000 to 2000 
integration time steps in order not to be overloaded with minutiae. 
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5. Subroutine RK4 

This subroutine was taken from the public domain collection of fortran 
subroutines held at the Naval Postgraduate School, It was written by Professor I. 
M. Ross and is a standard fourth order Runge-Kutta routine. It is used by 
passing the computed values for the variables and their derivatives to the routine 
via the variables X and XDOT. The routine actually requires four passes for 
the integration of one time step and during each of these four passes, the 
program goes back and recalculates the variables X and XDOT, The variable 
INDEX keeps track of the looping requirements for this routine. 

C. PROGRAM USAGE 

Proper operation of the program requires the manipulation of two input 
files. As previously mentioned, one file, AERO.DAT, contains the aerodynamic 
inputs and DATA.DAT contains the orbital dynamic inputs. The output file, 
OUT.DAT, contains the time history of ten important variables used in the 
analysis of the orbital plane changes. To set up a run, the user must first edit the 
AERO.DAT file to reflect the aerodynamics of the flight vehicle. The file 
AERO.DAT contains fifteen variables; each variable and the corresponding unit 
is explained in Table 1 . The second file DATA.DAT corresponds to the starting 
conditions of the flight vehicle. The file DATA.DAT contains eleven input 
variables all of which are contained in Table 2 along with the meanings of the 
variable and the necessary units. Once these two files have been set for the 
particular flight vehicle, the program ORBIT can be executed to obtain a time 
history output of the maneuver. 
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TABLE 1. INPUT FILE AERO.DAT 



Variable 


Meaning and Required Units 


Coefficient Cl 


Zero order coefficient for binomial 
curve fit of Cl vs alpha. 


Coefficient C2 


First order coefficient for binomial 
curve fit of Cl vs alpha. 


Coefficient C3 


Second order coefficient for binomial 
curve fit of Cl vs alpha. 


Coefficient C4 


Zero order coefficient for binomial 
curve fit of Cd vs alpha. 


Coefficient C5 


First order coefficient for binomial 
curve fit of Cd vs alpha. 


Coefficient C6 


Second order coefficient for binomial 
curve fit of Cd vs alpha. 


AOA for (CL/CD)max (AOPT) 


Angle of Attack for (CL/CD)max, 
Radians 


AOB for Bang Maneuver (AOB) 


Angle of Bang used during the 
Aerobang Maneuver, Radians. 


Density Exponent (N) 


Exponent of density within the heating 
rate model. Equation 4.5. 


Velocity Exponent (M) 


Exponent of velocity within the 
heating rate model. Equation 4.5. 


Reference Area (S) 


Flight vehicle reference area, used in 
lift and drag calculations, square meter 


Initial Mass (MASSO) 


Flight vehicle mass at beginning of 
maneuver, kg. 


Final Mass (FM) 


Final mass of vehicle at end of 
maneuver, kg. 


Specific Impulse (SPI) 


Specific Impulse of propulsion system, 
seconds. 


Thrust (N) 


Total thrust developed by propulsion 
system, or thrust used during 
Aerobang maneuver. Newtons. 
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TABLE 2. INPUT FILE DATA.DAT 



Variable 


Meaning and Required Units 


Case (CSE) 


Designate the type of model to be used. 


Radius (R) 


Geocentric radius of spacecraft body, 
units are in meters. 


Declination (0) 


Spherical coordinate angular 
measurement in radians, see Figure 12. 


Right Ascension ((j)) 


Spherical coordinate angular 
measurement in radians, see Figure 12. 


Speed (V) 


Speed of spacecraft, units are in meters 
per second. 


Flight Path Angle (y) 


Measure of spacecraft attitude 
referenced to local horizon, units in 
radians. 


Heading Angle (\j/) 


Measure of spacecraft attitude 
referenced to the earth's equator, units 
in radians. 


Begin Time (T) 


Start of simulation, units in seconds. 


End Time (TF) 


End of simulation, units in seconds. 


Time Interval 


Integration time step, units in seconds. 


Print Time Interval (TPI) 


Interval at which call is made to output 
subroutine, units in seconds. 



D. PROGRAM VALIDATION 

The program was validated in two ways, the first method involved testing the 
orbital mechanics portion of the program and the second method involved 
matching data output from a previously run problem. For the orbital mechanics 
validation all of the aerodynamic and thrust parameters were set to zero, while 
the control models were turned off and the flight velocity and altitude were set to 
a known orbital condition. The values output from the computer program were 
then compared to the calculated values. The comparison is summarized in Table 
3. As the data show, there is no appreciable difference between the calculated 
values and the computer output values. 
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TABLE 3. ORBITAL DYNAMIC VALIDATION 





Calculated Value | 


Program Value 


Apogee 


6675.0 km 


6675.0 km 


Apogee Velocity 


7.700 km/sec 


7.700 km/sec 


Period 


5370 sec 


5375 sec 


Eccentricity 


.007123 


.007129 


Perigee 


6580.5 kn 


6580.5 km 


Perigee Velocity 


7.810 km 


7.810 km 



The second phase of the validation involved matching the solution from the 



present computer program to a previously published solution using the same 
inputs. For this phase an aerocruise maneuver was simulated using the full 
capabilities of the program. 



TABLE 4. PROGRAM VALIDATION 





Value From Cervisi 
[Ref. 13] 


Program Value 


Velocity 


7.254 km/sec 


7.254 km/sec 


Radius of Orbit 


6447.5 km 


6447.5 km 


Initial/Final Mass 


4762.7/2948.4 kg 


4762.7/2948.2 kg 


Initial/Final AOB 


52.4/67.8 degrees 


52.37/67.8 degrees 


Inclination Change 


16.83 degrees 


16.87 degrees 



Reference 13 contains enough data to reconstruct the the necessary 



inputs for an accurate simulation. The data taken from the article was all in a 
non-dimensional format and it was necessary to take it back to a dimensional 
form for use in the program, ORBIT. Both outputs are summarized in Table 4 
and as can been seen, they match very closely. 
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VI. ANALYSIS 



Evaluation of the data gathered from modeling the atmospheric interaction 
of the two chosen flight vehicles will determine for what flight conditions the 
aerobang method is superior to the aerocruise. The analysis will demonstrate the 
similarity of the two control methods and the fact that the aerocruise model may 
be thought of as an aerobang maneuver constrained to a particular altitude and 
angle of attack. 

A. FORMULATION OF THE SOLUTION 

The objective was to conduct a numerical investigation into the efficiency of 
both the aerocruise and aerobang control laws. Efficiency of the maneuver is 
defined as the inclination change a maneuver produces for a given amount of fuel 
expended. The effects of three separate parameters on the efficiency of each 
maneuver were studied. The three parameters were the flight vehicle speed, the 
(L/D)niax of the vehicle, and the stagnation heating rate. In short, the maneuvers 
were conducted over a varying range of flight velocities for two different 
heating rates using first the ERV and then the MRRV. 

Flight vehicle velocities ranged from super- to subcircular velocities. The 
terms super- and subcircular are defined as faster or slower than the circular 
orbital velocity at that particular altitude. A parameter k can be defined as 

k = ^ (6.1) 

which is the ratio of the circular speed (Vc) to the flight speed. Thus, k>l is 
subcircular speed and k<l is supercircular speed. These two ranges are of 
interest because the control law for the aerobang maneuver dictates that at 
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supercircular speeds, the flight vehicle gains altitude, while at subcircular speeds, 
it might lose altitude. If the maneuver proves to be more efficient at 
supercircular speeds, the increase in altitude and velocity would only add to the 
efficiency by requiring less fuel to reorbit the spacecraft at termination of the 
atmospheric interaction. On the other hand, if the aerobang maneuver proves to 
be more efficient at subcircular speeds, a further analysis must be conducted to 
account for the loss in altitude and velocity. In this study, the altitude was 
maintained using two different methods. The first method involved giving the 
flight vehicle a small positive flight path angle, rather than a zero flight path 
angle; the upward tendency of the trajectory would offset the altitude loss by the 
maneuver. The second method was to tilt the angle of bank using some of the lift 
for the inclination change to keep the flight vehicle at a constant altitude. 

The maneuver simulations were conducted for two different heating-rate 
values in order to determine the effect different altitudes and velocities would 
have on the model of each maneuver. The heating rates were chosen to be 
0.9088* 10^ and 1.42*10^ Watts per square meter, referenced to a 1-ft radius 
sphere [Ref. 7: p. 490]. For the ERV, the maximum heating rate value would 
correspond a stagnation temperature of 4180 degrees Fahrenheit during a plane 
change maneuver lasting 1500 seconds [Ref. 14: pp. 278-9]. To effect an 
equitable comparison of the two flight vehicles using the same flight profiles, 
these heating rate values were also imposed on the MRRV, though the MRRV is 
designed to handle higher heating rates [Ref. 6: p. 516]. 

It is noted that the values of k selected for evaluation varied for each flight 
vehicle and for different heating rates. In order to evaluate the aerocruise and 
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aerobang maneuvers during similar flight condition the k values for the 
simulations are confined to a subset of flight speeds common to both maneuvers. 



B. SUPERCIRCULAR PROFILES 
1. High Heating Rate 
a. ERV 

A stagnation heating rate of 1.42*10^ watts per square meter was 
taken for the ERV. Simulations were run for two different supercircular speeds 
corresponding to values of k from .94 and .95. Each run was allowed to expend 
fuel up to 2% of the initial mass. The aerobang maneuver in either case did not 




Figure 17. Inclination Change Figure 18. Angle of Attack vs 
for ERV k=.94 Mass Fraction ERV k=.94 




Figure 19. Inclination Change Figure 20. Angle of Attack vs 
for ERV k=.95 Mass Fraction for ERV k=.95 



complete the full 2% mass expenditure. The maneuver was unable maintain the 
constant heating rate curve, thus terminating the flight profile. As shown in 
Figures 17 and 19 the aerobang maneuver did display a more efficient process of 
inclination change as compared to the aerocruise maneuver for supercircular 
speeds. 

For k=.94, the aerobang maneuver terminated at 1% mass fraction 
with a .226 degree inclination change, while the aerocruise achieved only .186 
degrees inclination change. At .75% mass fraction the inclination difference was 
at a maximum of .05 degrees. A strategy for maximizing inclination change 
would be to conduct multiple passes thought the atmosphere such that at each 
pass, the flight profile would be terminated at the point of maximum inclination 
change difference. The k=.95 flight profile displays similar characteristics to the 
k=.94 case except that the aerobang maneuver became less efficient at 
approximately .5% mass fraction of fuel. At lower values of fuel mass fraction 
the inclination change for the aerobang and aerocruise maneuvers were nearly 
indistinguishable from each other. 

No cases for k<.94 were run because of the inability of the 
aerocruise control law to maintain a constant heating rate while flying at the 
angle of attack corresponding to (L/D)niax< Th® aerocruise control law required 
angles of bank of 134 and 153 degrees for the k=.95 and the k=.94 cases 
respectively. For k=.93, the force required to balance the centrifugal force of 
the flight vehicle exceeded the lift generated, making it impossible for the flight 
vehicle to maintain the desired heating rate. In contrast, the aerobang maneuver 
could be executed under the restrictions/specifications of higher orbital speeds. 
A limitation encountered in the aerobang maneuver is that at faster speeds, the 
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angle of attack decreases to zero sooner, making the length of the maneuver 
shorter. 

Another difference in these two maneuvers is found in the time to 
complete the maneuver for a specified mass fraction of fuel. The aerobang 
maneuver is completed four times faster than the aerocruise maneuver for the 
same mass fraction burned. This observation has two implications: first, the 
aerobang maneuver can be completed closer to the node of the orbit giving the 
maneuver an additional increase in efficiency for long flight profiles, and 
second, the integrated heat load on the spacecraft is lower for the aerobang 
maneuver. The temperature of the vehicle will depend on the total heat brought 
into the vehicle minus the heat which can be dissipated away; if the heat load can 
be reduced, this would decrease the average temperature of the spacecraft. The 
reduced time of flight to complete an increased inclination change therefore 
implies that the aerobang maneuver could possibly allow flight profiles to be 
executed at a higher heating rate than the aerocruise maneuver, since the 
integrated heat load would be less, granting that the materials can handle the 
higher heating rate. 
b. MRRV 

The MRRV velocities were simulated from values of k from .99 to 
.98. No speeds greater than k=.98 were simulated due to the fact that the 
aerocruise maneuver could not generate the required lift at (L/D)max- All the 
simulations required fuel up to 2% of the spacecraft mass. Observing Figures 21 
through 23, it is clear that the aerobang maneuver is more efficient at producing 
an inclination change somewhere between k=.99 and k=.985. The thrust 
required to execute the maneuver using the aerobang method was 20 times 
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greater than the aerocruise method. This thrust differential is passed on to the 
time of flight required for each maneuver as the aerobang method required only 
19.3 seconds to expend the 2% fuel while it took over 400 seconds for each of 
the aerocruise cases. The thrust of the aerocruise maneuver is based on the drag 
encountered, while the MRRV has a high (L/D)max the actual drag force is low 
requiring only a small amount of thrust. 




Figure 21. Inclination Change Figure 22. Inclination Change 
for MRRV k=.99 for MRRV k=.985 




Figure 23. Inclination Change for MRRV k=.98 
2. Low Heating Rate 
a. ERV 

Simulations were conducted up to a mass fraction of 2% for fuel 
expended and a stagnation heating rate of 0.9088* 10^ watts per square meter 
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referenced to a 1 foot diameter sphere. This lower heating rate corresponds to a 
higher altitude and lower speed for the flight vehicle. The aerocruise simulation 
was unable to maintain altitude for k<.97 due to inadequate lift generated at 
(L/D)niax- The aerobang maneuver was able to expend all of the 2% mass 
fraction under these conditions. At k=.97, the aerobang maneuver once again 
proved to be more efficient at generating inclination change, with a total 
inclination change of .454 degrees at the end of the simulation. The point of 
maximum inclination difference between the aerocruise method and the aerobang 
method occurred at approximately 1% mass fraction, with the aerobang 
maneuver having a .265 degrees inclination change and the aerocruise maneuver 
having only .212 degrees inclination change. As an example, for k=.97 and a 
desired 10 degree inclination change, an aerobang maneuver would save 
approximately 490 kg of fuel using a multiple pass technique. Additionally, the 
time involved to conduct the aerocruise maneuver is almost nine times greater 
than the aerobang method. While there is a region of supercircular velocities 
where the aerocruise maneuver is more efficient than aerobang, it occurs nearer 
the circular velocity or k=l. 
b. MRRV 

The low heating rate flight profile for the MRRV was flown for 
speeds corresponding to k=.990 and k=.997. These velocities are very close to 
the circular orbital velocity. At the increased altitude and reduced velocity, as 
compared to the high heating rate simulation of section B.l.ii, the lift generated 
by the lifting body is reduced even further. Each of the aerocruise maneuvers 
require over 900 seconds to complete for a maximum inclination change of 1.1 
degrees at k=.997. The aerocruise maneuver exceeds the aerobang maneuver in 
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inclination change efficiency until sometime before the velocity exceeds the 
k=.995 level. Once again, the range in which the efficiency of the aerocruise 
method is superior is limited to a band close to the circular orbital velocity and 
beyond this band the aerobang method becomes increasingly more efficient. 

C. SUBCIRCULAR PROFILES 

Unlike the supercircular velocities the subcircular region requires more 
analysis to truly determine which method of inclination change is more efficient. 
Execution of the aerobang maneuver using a subcircular profile always appears 
to be more efficient but causes the spacecraft to lose altitude. This altitude loss 
must be made up for in order to effect an equitable comparison between the two 
maneuvers. Both methods of altitude maintenance discussed in the above section 
are utilized here to bring the vehicle back to nearly the same initial velocity and 
altitude. 

1. High Heating Rate 
a. ERV 

The ERV simulation was conducted for values of k=1.03 to k=l,l. 
Table 5 summarizes the data from these flight profiles. Column one of the table 
is the standard aerobang maneuver; column two is the aerobang maneuver using 
a positive flight path angle to return the flight vehicle to the original altitude at 
the end of the 2% mass fraction fuel expenditure; column three is the aerobang 
maneuver using an angle of bank to maintain a particular altitude, and column 
four is the aerocruise maneuver for which the comparison is made. Observing 
Table 5, the aerobang maneuver displays expected results: the efficiency of the 
maneuver remains higher than that of the aerocruise maneuver, but with an 
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TABLE 5. SUMMARY OF ERV FLIGHT PROFILES FOR HIGH 



HEATING RATE 





o 

II 


HEATING RATE 


1.42*10''6 WATTS/M ^^2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.046 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


83.4 deg 


73.3 


deg 


Inclination change 


0.93 deg 


0.754 deg 


0.765 deg 


0.849 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-124 m 


4 m 


-2 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


57.9 


sec 


fuel/degree inclination 


111.83 kg/deg 


137.93 kg/deg 


135.95 kg/deg 


122.5 


kg/deg 


% difference over 
Aerocruise 


8.71 % 


-12.60 % 


-10.98 % 


0 


% 




K= 1.05 


HEATING RATE 


1.42*10''6 WATTS/M''2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.072 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


80 deg 


65 


deg 


Inclination change 


1.05 deg 


0.779 deg 


0.794 deg 


0.821 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-196 m 


0 m 


-9 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


53.3 


sec 


fuel/degree inclination 


99.05 kg/deg 


133.5 kg/deg 


130.98 kg/deg 


126.67 


kg/deg 


% difference over 
Aerocruise 


21.80 % 


-5.39 % 


-3.40 % 


0 


% 




II 

o 


HEATING RATE 


1.42*10''6 WATTS/M''2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.097 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


77 deg 


58.3 


deg 


Inclination change 


1.166 deg 


0.794 deg 


0.823 deg 


0.787 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-264 m 


-2 m 


-15 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


49.1 


sec 


fuel/degree inclination 


89.19 kg/deg 


130.98 kg/deg 


126.37 kg/deg 


132.15 


kg/deg 


% difference over 
Aerocruise 


32.51 % 


0.89 % 


4.37 % 


0 


% 




K= 1.1 


HEATING RATE 


1.42*10-'6 WATTS/M''2 






AEROBANG 


AEROBANG^ 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.101 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


71 deg 


50.6 


deg 


Inclination change 


1.336 deg 


0.969 deg 


0.816 deg 


0.736 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-360 m 


-97 m 


1 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


43.6 


sec 


fuel/degree inclination 


77.84 kg/deg 


107.33 kg/deg 


127.45 kg/deg 


141.3 


kg/deg 


% difference over 
Aerocruise 


44.91 % 


24.04 % 


9.80 % 


0 


% 



attached altitude penalty. The altitude penalty increases as the percentage 
difference over the aerocruise maneuver increases. At k=1.03 the altitude loss is 
124 m while the difference over aerocruise is 8.7%; this increases at k=l.l to an 
altitude loss of 360 m with a difference of 44.9%. As shown in Figure 24 the 
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most effective method of altitude maintenance utilizes the angle of bank. The 
angle of bank method can be used over a much larger range of flight speeds than 
the flight path angle method; the flight path method cannot to be used beyond 
k=1.07 because no angle of attack solution exists for this combination of 
variables. From Figure 24, the aerobang maneuver that uses the angle of bank 
for altitude maintenance becomes more efficient than the aerocruise method at 
speeds less than k=1.06. As an example, for k=l.l a 20 degree inclination 
change using the aerobang method augmented with angle of bank would use 277 
kg of fuel less than aerocruise maneuver. 




Figure 24. Fuel Percent Difference over Aerocruise vs k Factor for 

ERV (High Heating Rate) 

Unlike the supercircular case in which the angle of attack during the 
aerobang maneuver steadily decreases to zero, in the subcircular case, it 
increases as the flight vehicle descends. The angle of attack could possibly 
exceed the stall angle of the flight vehicle, terminating the maneuver 
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prematurely. An unexpected advantage experienced during the angle-of-bank 
augmented aerobang maneuver is that the angle of attack remains constant 
throughout the maneuver, similar to the aerocruise maneuver, except the angle 
of attack is fixed at a higher value. This was observed in all cases where the 
aerobang maneuver was augmented with the angle of bank. Figures 25 and 26 
graphically show this concept. Figure 25 shows the inclination versus mass 
fraction for the aerobang maneuver at k=l.l and Figure 26 shows the inclination 
change for the same maneuver augmented with an angle of bank of 71 degrees. 





Figure 25. ERV k=l.l Angle of Figure 26. ERV k=l.l Angle of 
Attack Profile unaugmented Attack Profile AOB augmented 

b. MRRV 

The range of speeds at which the MRRV could be flown was 
significantly reduced when compared to the ERV. The parameter, k, ranged in 
value from 1.01 to 1.02; beyond 1.02 the speed was too slow to provide the 
necessary lift to maintain the aerocruise maneuver. This same effect was seen 
for the supercircular flight profiles. Table 6 summarizes the data obtained from 
the simulations conducted. Under these conditions, the efficiency of the flight 
path angle method of altitude maintenance was essentially identical to the angle of 
bank method. This may be attributed to the narrow range of flight velocities 
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sampled, all of which are closer to circular orbital velocity than the those 
sampled in the ERV case. 



TABLE 6. SUMMARY OF MRRV FLIGHT PROFILES FOR HIGH 



HEATING RATE 





o 

II 


HEATING RATE 


1.42*10''6 WATTS/M''2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.017 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


87.2 deg 


58.5 


deg 


Inclination change 


0.74 deg 


0.682 deg 


0.682 deg 


0.994 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-44 m 


1 m 


0 m 


0 


m 


Time required 


19.3 sec 


19.3 sec 


19.3 sec 


353.4 


sec 


fuel/degree inclination 


132.43 kg/deg 


143.70 kg/deg 


143.70 kg/deg 


98.59 


kg/deg 


% difference over 
Aerocruise 


-34.32 % 


-45.75 % 


-45.75 % 


0.00 


% 




K= 1.015 


HEATING RATE 


1.42*10''6 WATTS/M '‘2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.023 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


86.1 deg 


44.3 


deg 


Inclination change 


0.767 deg 


0.688 deg 


0.681 deg 


0.825 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-61 m 


0 m 


0 m 


0 


m 


Time required 


19.3 sec 


19.3 sec 


19.3 sec 


346 


sec 


fuel/degree inclination 


127.77 kg/deg 


142.44 kg/deg 


143.91 kg/deg 


118.79 


kg/deg 


% difference over 
Aerocruise 


-7.56 % 


-19.91 % 


-21.15 % 


0 


% 




K= 1.02 


HEATING RATE 


1.42*10''6 WATTS/M '‘2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.029 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


85 deg 


26.4 


deg 


Inclination change 


0.794 deg 


0.694 deg 


0.69 deg 


0.542 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-78 m 


-1 m 


1 m 


0 


m 


Time required 


19.3 sec 


19.3 sec 


19.3 sec 


338.8 


sec 


fuel/degree inclination 


123.43 kg/deg 


141.21 kg/deg 


142.03 kg/deg 


180.81 


kg/deg 


% difference over 
Aerocruise 


31.74 % 


21.90 % 


21.45 % 


0 


% 



Figure 27 shows the break even point for the efficiency comparison 
of aerobang to aerocruise to be at approximately k=1.0175; slower than this 
speed the aerobang method is more efficient in providing an inclination change. 
At k=1.02, the difference in fuel of the aerobang method over the aerocruise 
method exceeds 20 percent. This high difference is attributed to the excessive 
angle of bank required by the aerocruise maneuver to maintain the flight profile. 
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The aerocruise maneuver for this flight condition is at an angle of bank of 26 
degrees, allowing a large portion of the lift to be directed away from the orbit 
normal and unable to effect the inclination change. 




Figure 27. Fuel Percentage Difference over Aerocruise vs k Factor 

for MRRV (High Heating Rate) 

2. Low Heating Rate 
a. ERV 

The range of flight velocities over which the aerocruise maneuver 
could operate was narrower from those of the high heating rate. Speeds slower 
than k=1.05 required an unattainable angle of attack for the aerocruise 
maneuver. Table 7 summarizes the data gathered for the given set of flight 
conditions. Comparing the values of specific inclination change (ie, fuel/degree 
inclination) the two aerobang augmented methods were nearly identical. A 



54 



TABLE 7. SUMMARY OF ERV FLIGHT PROFILES FOR LOW 



HEATING RATE 





K= 1.02 


HEATING RATE .9088*10''6 


WATTS/M^2 




AERCffiANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.057 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


84.5 deg 


65.3 


deg 


Inclination change 


1.675 deg 


1.303 deg 


1.3 deg 


1.603 


deg 


Mass Fraction used 


4% 


4% 


4% 


4% 




Change in Altitude 


-358 m 


-6 m 


3 m 


0 


m 


Time required 


40.9 sec 


40.9 sec 


40.9 sec 


239.1 


sec 


fuel/degree inclination 


124.2 kg/deg 


159.6 kg/deg 


160 kg/deg 


129.8 


kg/deg 


% difference over 


4.31 % 


-22.96 % 


-23.27 % 


0.00 


% 


Aerocruise 














K= 1.03 


HEATING RATE 


.9088*10*6 WATTS/M*2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.044 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


82.2 deg 


45.3 


deg 


Inclination change 


0.788 deg 


0.655 deg 


0.653 deg 


0.635 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-124 m 


0 m 


-1 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


140.9 


sec 


fuel/degree inclination 


131.98 kg/deg 


158.78 kg/deg 


159.26 kg/deg 


163.78 


kg/deg 


% difference over 


19.42 % 


3.05 % 


2.76 % 


0 


% 


Aerocruise 














K= 1.04 


HEATING RATE 


.9088*10*6 WATTS/M*2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.057 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


79.8 deg 


28.9 


deg 


Inclination change 


0.836 deg 


0.665 deg 


0.658 deg 


0.446 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-160 m 


-2 m 


1 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


135.2 


sec 


fuel/degree inclination 


124.40 kg/deg 


156.39 kg/deg 


158.05 kg/deg 


233.18 


kg/deg 


% difference over 


46.65 % 


32.93 % 


32.22 % 


0 


% 


Aerocruise 














K= 1.05 


HEATING RATE 


.9088*10*6 WATTS/M*2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.069 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


78 deg 


3.7 


deg 


Inclination change 


0.873 deg 


0.668 deg 


0.661 deg 


0.135 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-189 m 


0 m 


1 m 


0 


m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


130.8 


sec 


fuel/degree inclination 


119.13 kgid&g 


155.69 kg/deg 


157.34 kg/deg 


770.37 


kg/deg 


% difference over 


84.54 % 


79.79 % 


79.58 % 


0 


% 


Aerocruise 













demonstrated disadvantage of the flight path angle method of augmentation is 
that the angle of attack for each simulation exceeded 40 degrees, beyond the 
flight vehicle’s operating envelope. This excessive angle of attack did not occur 
for any of the angle of bank augmented simulations and while the angle of attack 
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was high it did not exceed the threshold of the model. Additionally, from Figure 
28 it is observed that the aerobang maneuver is more efficient at k=1.03 or 
greater. Once again the aerobang maneuver has a higher specific inclination 
change compared to the aerocruise maneuver as the velocity differs more from 
the circular orbital speed. The time for the two methods to complete the flight 
profile differed greatly: the aerobang method once again was much faster than 
the aerocruise method, completing the maneuver approximately 6.5 times faster. 
For short duration atmospheric interactions, a time differential of this nature 
may not play an important roll, however, if the atmospheric interactions were of 




Figure 28. Fuel Percentage Difference over Aerocruise vs k Factor 

for MRRV (Low Heating Rate) 
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longer duration, the aerocruise maneuver may approach or exceed a quarter 
period of its orbit, in which case the inclination change caused by any force in 
the orbit normal direction is reduced. 
b. MRRV 

The last set of simulations discussed will be that of the MRRV at a 
low heating rate value. The MRRV range of flight speeds are very narrow at 
this heating rate: from k= 1.002 to k= 1.007. The aerocruise maneuver cannot be 
performed at speeds slower than this. As previously stated, the narrow range of 



TABLE 8. SUMMARY OF MRRV FLIGHT PROFILES FOR LOW 



HEATING RATE 





K= 1.002 


HEATING RATE 


.9088*10^6 WATTS/M^2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.006 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


89 deg 


62.5 


deg 


Inclination change 


0.637 deg 


0.619 deg 


0.621 deg 


0.984 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-15 m 


0 m 


-1 m 


0 


m 


Time required 


19.3 sec 


19.3 sec 


19.3 sec 


891 


sec 


fuel/degree inclination 


1 . 09 kg/deg 


158.32 kg/deg 


157.81 kg/deg 


99.59 


kg/deg 


% difference over 
Aerocruise 


80.62 % 


-58.97 % 


-58.45 % 


0.00 


% 




K= 1.005 


HEATING RATE 


.9088*10*6 WATTS/M *2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.009 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


88.2 deg 


39.3 


deg 


Inclination change 


0.651 deg 


0.625 deg 


0.623 deg 


0.714 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-26 m 


-1 m 


0 m 


0 


m 


Time required 


19.3 sec 


19.3 sec 


19.3 sec 


879.6 


sec 


fuel/degree inclination 


1.09 kg/deg 


156.80 kg/deg 


157.30 kg/deg 


137.25 


kg/deg 


% difference over 
Aerocruise 


85.94 % 


-14.24 % 


-14.61 % 


0 


% 




K= 1.007 


HEATING RATE 


.9088*10*6 WATTS/M*2 






AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.013 deg 


0 deg 


0 


deg 


Initial AOB 


90 deg 


90 deg 


67.7 deg 


12.8 


deg 


Inclination change 


0.661 deg 


0.623 deg 


0.623 deg 


0.29 


deg 


Mass Fraction used 


2% 


2% 


2% 


2% 




Change in Altitude 


-34 m 


0 m 


0 m 


0 


m 


Time required 


19.3 sec 


19.3 sec 


19.3 sec 


872.2 


sec 


fuel/degree inclination 


1.09 kg/deg 


157.30 kg/deg 


157.30 kg/deg 


337.93 


kg/deg 


% difference over 
Aerocruise 


94.29 % 


53.45 % 


53.45 % 


0 


% 
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speeds is attributed to the fact that the MRRV was designed to operate 
aerodynamically at lower altitudes and higher heating rates than the value used. 
For comparison to the ERV, the range is adequate. Table 8 is the summary of 
each set of data and Figure 29 graphically represents the data contained in the 
table. From Figure 29 the aerobang method of inclination change becomes more 
efficient at speeds slower than k= 1.0055. The aerocruise maneuver required 
nearly 900 seconds to be completed as compared to the 19 seconds for the 
aerobang maneuver. At 900 seconds the aerocruise maneuver approaches a 




Figure 29. Fuel Percentage Difference over Aerocruise vs k Factor 

for MRRV (Low Heating Rate) 
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quarter of an orbital period, beyond the 900 seconds shown, the inclination 
change per kilogram of fuel falls off to become zero at the apex of the orbit. 
For this maneuver, the flight vehicle would need to roll over to bring the lift 
vector to the other side of the orbit or be subject to a decreasing inclination. 

D. MASS FRACTION EFFECTS 

All of the above analysis was done for fuel expenditure of 2% of spacecraft 
mass fraction. Because larger inclination changes may be desired, an analysis 
was done on the possible advantages of expending more than 2 % mass fraction of 
fuel. The effect of longer atmospheric interactions may cause the maneuver to 
approach the apex of the orbit and decrease the inclination change. The 
aerobang method is faster at completing the desired inclination change as 
compared to the aerocruise method. A conclusion can be drawn that the 
aerobang method should become increasingly better than the aerocruise method 
as the mass fraction increases. Table 9 shows the effect of increasing mass 
fraction for the ERV with a flight profile of k=1.03 and a high heating rate. 
Both methods of altitude maintenance were used and compared to the aerocruise 
method. The flight path angle method was terminated beyond a mass fraction of 
4% because the vehicle could not be flown at the constant heating rate. The 
angle of bank method of altitude maintenance worked throughout the full range 
of mass fraction. Figure 30 shows the percentage difference over aerocruise vs 
the mass fraction. The graph indicates that initially the mass fraction has little 
effect on the percentage difference over the aerocruise model. As more mass 
fraction is burned, the effect was as expected; the aerocruise method lost some of 
its advantage because it approaches the apex of the orbit. For this particular 
flight profile, the aerocruise maneuver took approximately three times longer 
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than the aerobang maneuver. For flight profiles which have a higher ratio of 
maneuver times, the mass fraction advantage effect would occur at lower mass 
fractions than for the case presented. 



TABLE 9. EFFECT OF MASS FRACTION ON INCLINATION 



CHANGE EFFICIENCY 





K= 1.03 


HEATING RATE 


1.42*10''6 WATTS/M*2 




AEROBANG 


AEROBANG 


AERC©ANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.046 deg 


0 deg 


0 deg 


Initial AOB 


90 deg 


90 deg 


83.4 deg 


73,3 deg 


Inclination change 


0.93 deg 


0.754 deg 


0.765 deg 


0.849 deg 


Mass Fraction used 


2% 


2% 


2% 


2% 


Change in Altitude 


.124 m 


4 m 


-2 m 


0 m 


Time required 


20.5 sec 


20.5 sec 


20.5 sec 


57.9 sec 


fuel/degree inclination 


111.83 kg/deg 


137.93 kg/deg 


135.95 kg/deg 


122.5 kg/deg 


% difference over 


8.71 % 


-12.60 % 


-10.98 % 


0 % 


Aerocruise 












o 

II 


HEATING RATE 


1.42*10''6 WATTS/M*2 




AEROBANG 


AEROBANG 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.08 deg 


0 deg 


0 deg 


Initial AOB 


90 deg 


90 deg 


83.4 deg 


73.3 deg 


Inclination change 


2.193 deg 


1.547 deg 


1.546 deg 


1.716 deg 


Mass Fraction used 


4% 


4% 


4% 


4% 


Change in Altitude 


-516 m 


-29 m 


-6 m 


0 m 


Time required 


40.9 sec 


40.9 sec 


40.9 sec 


115.7 sec 


fuel/degree inclination 


94.8 kg/deg 


134.45 kg/deg 


134.36 kg/deg 


121.2 kg/deg 


% difference over 


21,78 % 


-10.93 % 


-10.86 % 


0 % 


Aerocruise 












o 

II 


HEATING RATE 


1.42*10''6 WATTS/M*2 




AEROBANG 


AEROBANG* 


AEROBANG 


AEROCRUISE 


Initial Gamma 


0 deg 


0.097 deg 


0 deg 


0 deg 


Initial AOB 


90 deg 


90 deg 


83.42 deg 


73.3 deg 


Inclination change 


8.684 deg 


5.964 deg 


3.945 deg 


4.412 deg 


Mass Fraction used 


10% 


10% 


10% 


10% 


Change in Altitude 


.4092 m 


-1760 m 


4 m 


0 m 


Time required 


102 sec 


102 sec 


102 sec 


288.5 sec 


fuel/degree inclination 


59.65 kg/deg 


86.85 kg/deg 


131.3 kg/deg 


117.5 kg/deg 


% difference over 


49.23 % 


26.09 % 


-11.74 % 


0 % 


Aerocruise 












K= 1.03 


HEATING RATE 


1.42*10*6 WATTS/M *2 




AEROBANG* 


AEROBANG* 


AEROBANG 


AEROCRUISE 


Initial Gamma 


deg 


deg 


0 deg 


0 deg 


Initial AOB 


deg 


deg 


83.77 deg 


73.3 deg 


Inclination change 


deg 


deg 


19.29 deg 


20.095 deg 


Mass Fraction used 






40% 


40% 


Change in Altitude 


m 


m 


0 m 


0 m 


Time required 


sec 


sec 


408.2 sec 


1153.8 sec 


fuel/degree inclination 


kg/deg 


kg/deg 


107.5 kg/deg 


103.2 kg/deg 


% difference over 


% 


% 


-4.17 % 


0 % 


Aerocruise 
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Figure 30. Fuel Percentage Difference over Aerocruise vs Mass 
Fraction Fuel Expended for ERV (k=1.03) 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



1. CONCLUSIONS 

As demonstrated in Chapter VI, there exists a small range of flight speeds in 
which the aerocruise maneuver exhibits a greater inclination change for a given 
amount of fuel than does the aerobang maneuver. This range of flight speeds is 
always centered around circular orbital speed. Figures 31 and 32 graphically 
shows the flight range where (within the hash marks), the aerocruise maneuver is 
predicted to be more efficient. However, use of the aerobang maneuver is 
predicted to be more efficient as the speeds become either more super- or more 
subcircular (regions beyond the hash marks). These regions of superior 
aerocruise efficiency, while bracketing the circular orbit speed are not 
symmetric in nature and tend to be narrower in the supercircular flight regime. 




k factor 



Figure 31. Efficiency Regions for the ERV 
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Figure 32. Efficiency Regions for the MRRV 
As discussed in Chapter IV, Section A, the aerocruise maneuver was unable to 
maintain the desired profile at speeds just possible the region of superior 
efficiency. This result was expected and shows that the the aerobang maneuver 
has a greater range of operating flight speeds as compared to the aerocruise 
maneuver. 

A higher maneuver efficiency exists for a vehicle with a greater maximum 
lift-to-drag ratio while executing an aerobang maneuver. The data from the 
ERV and the MRRV were compared for similar flight profiles, for the case 
k=1.03 and Q= 1.42*10^6 w/m^2. In order to analyze the effect of the lift-to- 
drag ratio on the inclination change, the effects from factors not associated with 
the lift-to-drag ratio must be eliminated. A major component which directly 
influences the amount of inclination change achieved is the lifting force 
generated by the flight vehicle during the maneuver. Lift generated by the flight 
vehicle is the product of the dynamic pressure on the vehicle, the coefficient of 
lift, and the reference wing area. By operating the two vehicles at the same 
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heating rate and flight speed, the dynamic pressure remains the same for the two 
simulations. The reference wing area of the ERV and the MRRV differ and 
must be taken into account in order to isolate the effect the lift-to-drag ratio has 
on the inclination change. 

A scaling factor equivalent to the ratio of the wing areas of ERV and the 
MRRV is introduced. This scaling factor is used to account for the inclination 
change lost due to the lower wing area of the MRRV. From Chapter III, the 
wing area of the MRRV is 28% lower than that of the ERV, yielding a scaling 
factor of 1 .4. Prior to scaling, the inclination change for the MRRV under these 
flight conditions is 0.693 degrees for a 2% fuel consumption. The ERV attained 
an inclination change of 0.765 degrees under identical flight conditions. Scaling 
the MRRV inclination change with the 1.4 factor yields an inclination change of 
0.970 degrees, clearly higher than that for the ERV. 

As can be seen in Figure 32 the aerocruise-superior flight speed range of 
operations for the MRRV dramatically narrowed; this may be attributed to the 
the decreased lift force generated by the smaller coefficient of lift rather than the 
higher lift-to-drag ratio of the MRRV. The time of flight for the aerobang 
maneuver was unaffected by the change in the lift-to-drag ratio while for the 
aerocruise maneuver, the time required to bum a 2% mass fraction of fuel 
increased for the MRRV. This result was not unexpected because the reduced 
drag coefficient of the higher lift-to-drag ratio vehicle is reflected in the thrust 
required to maintain altitude. 

Operating at a subcircular speed, the aerobang maneuver would normally 
results in a loss of altitude. In order to compare this maneuver with the 
aerocruise maneuver, which loses no altitude, two methods of altitude control 
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were used (Section VI.A). Of the two methods of altitude maintenance, 
augmentation of the angle of bank proved to be superior to modification of the 
initial flight path angle. While there was generally little difference in efficiency 
from using either method. Figure 24 clearly shows the angle of bank method to 
be more efficient than the flight path angle method. Additionally, the useful 
range of operation of the flight path angle method is much smaller than that of 
the angle of bank method. In all subcircular cases, there existed speeds too low 
for the flight path angle method to maintain a flight profile at a constant heating 
rate. 

From a fuel efficiency point of view, the percentage difference in fuel 
burned per unit angle of inclination change between the aerobang and aerocruise 
methods remains unaltered until the time required to complete the maneuver 
approaches a quarter of an orbital period for the aerocruise maneuver. This 
means that a low drag environment (drag force « maximum thrust), the 
aerobang maneuver will perform better than the aerocruise as more fuel is 
expended. The reason for this is that the thrust for the aerobang maneuver is not 
dependent on the drag environment, while the thrust for the aerocruise maneuver 
is modulated to counteract the drag. 

2. FOLLOW ON TOPICS 

Simply put, during the maneuver, the aerobang method modulates the angle 
of attack and the aerocruise method modulates the angle of bank. As shown 
within this thesis, there exist regions where one method claims superiority over 
the other. Following this line of thinking, a method of control which modulates 
both the angle of attack and the angle of bank to realize the best of each 
maneuver may be of greater efficiency than either of the two methods studied 
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here. This would require more sophisticated control laws than those discussed 
above. 

Another topic of interest would be to determine whether to use the angle of 
bank method of maintaining altitude during an aerobang maneuver, or to accept 
the energy loss normally associated with the subcircular aerobang maneuver and 
reestablish the orbit with an additional bum after exiting the atmosphere. It may 
take less fuel to accept the altitude loss from a true aerobang maneuver than to 
augment the maneuver in order to maintain a specified altitude. 

Comparison of the aerocruise data for two flight vehicles operating at 
similar flight conditions produced a curious result. It is well established that a 
vehicle with the higher lift-to-drag ratio produces the more efficient maneuver, 
this was not observed for the aerocruise maneuver and further investigation is 
warranted. This may perhaps be attributed to the fact that although the 
aerocruise maneuver operates at the maximum lift-to-drag ratio for the vehicle, 
the actual values of the drag and lift may be low. If the lift force is low, then the 
plane-change achieved will also be low for a given arc length regardless of the 
maximum lift-to-drag ratio. 

Finally, there is reason to investigate the unsymmetric nature of the 
aerocruise- superior region of operations (Figure 32). For both flight vehicles, 
the supercircular region of the aerocruise-superior area is narrower than the 
subcircular region. The flight regime where the aerocruise maneuver is 
superior was expected to center around the circular orbital speed, but in a 
symmetric nature. The reason why this occurred is unclear at this time and 
should be addressed in a follow-on study. 
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APPENDIX A 



C 



C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



PROGRAM ORBIT 
IMPLICIT REAL*8(A-H,L-Z) 

DIMENSION X(6),XDOT(6),C(6) 

CHARACTER* 10 CSE 

PROGRAMMER : THOMAS P. SPRIESTERBACH 

DATE : 15 AUG 91 

SYSTEM : VAX/FORTRAN 



VARIABLE LIST 



I 

INDEX 

J 

KOUNT 

DUM 

OLDCO 

CHANGE 

NE 

CSE 

T 

TF 

H 

TPI ' 
X(l) 

X(2) 

X(3) 

X(4) 

X(5) 

X(6) 

XDOTO 

C(l) 

C(2) 

C(3) 

C(4) 

C(5) 

C(6) 

BETA 

GO 

G 

HO 

RHOO 

MU 

ALPHA 



COUNTING INDEX 

COUNTING INDEX FOR RUNGE-KUTTA ROUTINE 
COUNTING INDEX 

COUNTING INDEX FOR OUTPUT DETERMINATION 
DUMMY VARIABLE FOR CONVENIENCE 
CONTROL VARIABLE FOR CONVERGENCE ON AOA 
INCREMENTAL CHANGE IN AOA USING NEWTONS METHOD 
NUMBER OF EQUATIONS OF MOTION 
CASE EITHER AEROCRUISE OR AEROBANG 
TIME (SEC) 

FINAL TIME (SEC) 

INCREMENT OF TIME [RUNGE-KUTTA ROUTINE](SEC) 

TIME PRINT INTERVAL (SEC) 

RADIUS (METERS) 

THETA (RADIANS) 

PHI (RADIANS) 

VELOCITY (M/S) 

FLIGHT PATH ANGLE GAMMA (RADIANS) 

PSI (RADIANS) 

DERIVATIVES OF THE ABOVE SEX VARIABLES 
ZEROTH ORDER COEFnCEENT FOR CL EQUATION 
FIRST ORDER COEFFICIENT FOR CL EQUATION 
SECOND ORDER COEFFICIENT FOR CL EQUATION 
ZEROTH ORDER COEFHCIENT FOR CD EQUATION 
FIRST ORDER COEFFICIENT FOR CD EQUATION 
SECOND ORDER COEFHCIENT FOR CD EQUATION 



DENSITY MODEL EXPONENT 
GRAVITY AT EARTH SURFACE 
LOCAL GRAVITATIONAL ATTRACTION 
REFERENCE ALTITUDE FOR DENSITY MODEL 
REFERENCE DENSITY FOR DENSITY MODEL 
GRAVITATIONAL CONSTANT 
ANGLE OF ATTACK 



(METERS) 
(9.806 M'^2/S) 
(M'^2/S) 
(METERS) 
(KG/M'^S) 

(RADIANS) 
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c 


ALFA 


ANGLE OF ATTACK 


(DEG) 


c 


AOPT 


ANGLE OF ATTACK AT MAX OVCD 


(RADIANS) 


c 


MASS 


MASS OF SPACECRAFT 


(KG) 


c 


MASSO 


INITIAL MASS OF SPACECRAFT 


(KG) 


c 


SPI 


SPECinC IMPULSE OF POWER PLANT 


(SEC) 


c 


R 


RADIUS 


(M) 


c 


V 


VELOCITY 


(M/S) 


c 


RADIUS 


RADIUS 


(KM) 


c 


VELOCITY 


VELOCITY 


(KM/S) 


c 


GAMMA 


FLIGHT PATH ANGLE 


(DEG) 


c 


ANGM 


ANGULAR MOMENTUM 


(M'^2/S) 


c 


SIGMA 


ANGLE OF BANK 


(DEG) 


c 


OINC 


INCLINATION 


(DEG) 


c 


DELTA 


INTERMEDIATE ANGLE 


(DEG) 


c 


OMEGA 


LONGITUDE OF THE ASCENDING NODE 


(DEG) 


c 


AS 


TANGENTIAL ACCELERATION 


(M/S'^2) 


c 


AR 


RADIAL ACCELERATION 


(M/S'^2) 


c 


AW 


BINORMAL ACCELERATION 


(M/S'^2) 


c 


THR 


THRUST 


(N) 


c 


AOB 


ANGLE OF BANK 


(RADIANS) 


c 


RHO 


DENSITY 


(KG/M^3) 


c 


S 


REFERENCE AREA OF SPACECRAFT 


(M'^2) 


c 


N 


DENSITY COEFFICIENT FOR HEATING RATE EQUATION 


c 


M 


VELOCITY COEFFICIENT FOR HEATING RATE EQUATION 


c 


HEAT 


STAGNATION POINT HEATING RATE 


(WATTS/M^2) 


c 


QS 


DYNAMIC PRESSURE TIMES THE REF AREA 


(N*M'^2) 


c 


CD 


COEFFICIENT OF DRAG 




c 


CL 


COEFHCIENT OF LIFT 




c 


LIFT 


LIFT ON SPACECRAFT 


(N) 


c 


DRAG 


DRAG ON SPACECRAFT 


(N) 



c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C DEFINED FUNCTION FOR ACCEL DUE TO GRAVITY 
C 

G(R)=MU/(R*R) 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C OPEN TWO INPUT FILES DATA.DAT - ORBITAL FARMS 

C AERO.DAT - VEHICLE FARMS 

C ONE INPUT FILE OUT.DAT - TABLE OF VALUES 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

OPEN(10,nLE=’DATA’,STATUS=’OLD’) 

OPEN(l 2, FILE=’ AERO', STATUS ='OLD') 

OPEN( 1 3,FILE=’OUr,STATUS=’NEW') 

READ( 1 0, 1 )CSE,(X(I),I= 1 ,6),T,TF,H,TPI 
1 FORM AT(/,/,20X, A 1 0, 1 0(/,20X,D 1 3.7)) 
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2 



RE AD( 1 2,2)(C(I),I=1 ,6), AOPT, AOB ,N,M,S ,MAS SO,FM,SPI,THR 
FORMAT(/,/yy, 1 5(/,20X,D 13.7)) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C PRINT CASE TO SCREEN FOR VISUAL VERIFICATION OF SELECTION 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

WRITE(6,*)CSE,'*' 

WRITE(6,*) 

c 

C NUMBER OF EQUATIONS TO INTEGRATE 
C X(1)....RADIUS 

C X(2)....THETA (SPHERICAL COORD. PARAMETER) 

C X(3)....PHI (SPHERICAL COORD. PARAMETER) 

C X(4)....VELOCITY 
C X(5)....FLIGHT PATH ANGLE GAMMA 
C X(6)....PSI 
C 

C INITIALIZATION OF CONSTANTS 
C 

NE=6 

PI=3.14159265359D+0 

MU=3.986012D+14 

G0=9.806D+0 

MASS=MASSO 

C 

C MAIN PROGRAM 
C 

INDEX=0 

KOUNT=l 

C 

C INITIAL CALL TO GET ACCELERATIONS AND CONTROL VARIABLES 
C FOR HRST OUTPUT 
C 

CALL CNTRL(DRAG,LIFT,THR, ALPHA, AOB,X, MASS, S,C,N,M, 
♦RHO,CSE,AOPT) 

CALL ACEL(AS,AR,AWX>RAG,LIFT,THR,ALPHA,AOB,MASS) 

CALLHDR(X(1),X(4),MASS,THR,X(5),CSE) 

CALLWRT(X,T,THR,MASS,ALPHA,AOB,RHO,N,M) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C THIS IS THE MAIN BLOCK OF THE PROGRAM CALLING RRST THE 
C CONTROL SUBROUTINE TO DETERMINE THE CONTROL METHOD AND 
C OUTPUT VARIABLES. NEXT ON TO THE ACCELERATION ROUTINE TO 
C CALCULATE THE RELATIVE ACCELERATION TO THE VEHICLE, THEN ON 
C TO THE DIFFERENTIAL EQUATIONS OF MOTION AND FINALLY A 
C FOURTH ORDER RUNGE-KUTTA ROUTINE. THIS BLOCK IS COMPUTED 
C AT EACH TIME INCREMENT H. 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

100 CALL CNTRL(DRAG, LIFT, THR,ALPHA,AOB,X, MASS, S,C,N,M, 
*RHO,CSE,AOPT) 

CALL ACEL(AS,AR,AWJ)RAG,LIFT,THR,ALPHA,AOB,MASS) 

CALL ORB(X,XDOT,AS,AR,AW,G(X(l))) 

CALL RK4(T,X,XDOT,NE,H,INDEX) 

IFONDEX .NE. 0) GO TO 100 

C DETERMINE IF TIME TO PRINT TO OUTPUT USING TIME PRINT 
C INTERVAL 

200 IF(KOUNT .LT. IDNINT(TPI/H)) GO TO 300 

CALLWRT(X,T,THR,MASS,ALPHA,AOB,RHO,N,M) 

C PRINT STATUS TO TERMINAL 

WRrrE(6,7)T,ALPHA,AOB 
7 FORMAT('+’,F8.1,lX,F6.3,lX,F6.3) 

KOUNT=0 

300 KOUNT=KOUNT+l 
C 

C CHANGE MASS OF S/C 
C 

MASS=MASS-THR*H/(SPI*G0) 

C 

C TWO TESTS EITHER CAN STOP THE PROGRAM. TEST FOR MASS LESS 
C THAN FINAL MASS AND ALSO IF TIME IS GREATER THAN FINAL TIME 
C 

IF(MASS .GT. FM .AND. T .LE. TF) GO TO 100 
C 

C PRINT OUTPUT AGAIN IF MASS SWITCH ENDED PROGRAM 
C 

IF(MASS .LE. FM) CALL WRTCX.T.THR, MASS, ALPHA, AOB,RHO,N,M) 

C 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C THIS SUBROUTINE CALCULATES THE CONTROL OF THE FREE 
C VARIABLES DURING THE ATMOSPHERIC PENETRATION 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE CNTRL(DRAG,LIFT,THR,ALPHA,AOB,X,MASS,S,C,N,M 
*,RHO,CSE,AOPT) 

IMPLICIT REAL*8(A-HJ.-Z) 

CHARACTER*10 CSE 
DIMENSION X(6),C(6) 

C 

MU=3.986012D+14 

GAMMA=X(5) 
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R=X(1) 

V=X(4) 



COEFHCIENTS FOR AN EXPONENTIAL DENSITY MODEL 
REFERENCE J MEASE 1976 US STANDARD ATMOSPHERE 
RANGE 50KM TO 120KM 

BETA=1.41D-4 

RHO0=3.0968D-4 

R0=6.438D+6 

G=MU/(R*R) 

RHO=RHOO*EXP(-BETA*(R-RO)) 

QS=.5D-K)*RHO*V*V*S 



THE FIRST CASE IS FOR THE AEROB ANG CONTROL LAW AND THE 
SECOND IS FOR THE AEROCRUISE CONTROL LAW. FOR AEROB ANG THE 
AOB IS SET AND THE ANGLE OF ATTACK IS CONTROLLED TO FLY AT A 
CONSTANT HEATING RATE, THE ALTITUDE AND VELOCITY ARE 
ALLOWED TO FLOAT 

IF (CSE.EQ.'AEROBANG ') THEN 

GUESS ALPHA INITIALLY FOR NEWTON APPROXIMATION ROUTINE 

ALPHA=(C(5)*QS+SQRT(C(5)**2+4.*QS*C(6)*(THR-QS*C(4)))) 

*/(2.*C(6)*QS) 

USE NEWTON APPROXIMATION METHOD TO CONVERGE ON ALPHA FOR 
AEROBANG USE A WHILE STRUCTURE FOR THE CONVERGENCE OF 
ALPHA SET OLDCO AND CHANGE EQUAL TO THE VALUES BELOW TO 
ENSURE THAT THE PROGRAM SWITCHES ARE NOT INITIALLY TRIPPED 

OLDCO=1X)+10 

CHANGE=L 

DO WHILE (ABS(CHANGE) .GT. l.D-3) 

CL=C(1)+C(2)*ALPHA+C(3)*ALPHA**2 

CD=C(4)+C(5)*ALPHA+C(6)*ALPHA**2 

LIFT=CL*QS 

DRAG=CD*QS 

THESE THREE EQUATIONS REPRESENT THE FUNCTION AND ITS 
DERIVATIVE WITH RESPECT TO ALPHA FOR NEWTON’S METHOD OF 
APPROXIMATION 

COEFF=(THR*COS(ALPHA)-DRAG)- 

*MASS*SIN(GAMMA)*(G+BETA*(N/M)*V*V) 
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COEFP=-THR*SIN(ALPHA)-(C(5)+2*C(6)*ABS(ALPHA))*QS 

C 

CHANGE=COEFF/COEFP 

ALPHA=ALPHA-CHANGE 
IF (ABS(OLDCO)-ABS(CHANGE).LE.O.) STOP 
OLDCO=CHANGE 
END DO 
C 

ELSE 

C 

C THIS IS THE CONTROL ROUTINE FOR AEROCRUISE. THRUST COUNTER 
C BALANCES THE DRAG AND THE ANGLE OF BANK IS SET TO COL i ROL 
C THE ORBIT AT CONSTANT ALTITUDE 
C 

ALPHA=AOPT 

C 

CL=C( 1 )+C(2)* ALPHA+C(3)* ALPHA**2 
CD=C(4)+C(5)*ALPHA+C(6)*ALPHA**2 
C 

LIFT=CL*QS 

DRAG=CD*QS 

C 

THR=DRAG/COS(ALPHA) 

C FORTRAN STOP IF IT IS IMPOSSIBLE TO FLY AN AEROCRUISE PROFILE 
C THIS OCCURS WHEN THE PROGRAM TRIES TO FIND AN AOB IN WHICH 
C THE COSINE OF ms GREATER THAN ONE. 



IF(ABS(MASS*(((G*R-V**2)/R)/(THR*SIN(ALPHA)+LIFT))).GT.l.)STOP 

AOB=ACOS(MASS*(((G*R-V**2)/R)/(THR*SIN(ALPHA)+LIFT))) 

C 

ENDIF 

C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C THIS SUBROUTINE DEFINES THE ACCELERATIONS TANGENTIAL, 

C PERPENDICULAR AND RADIAL TO THE ORBITAL PLANE 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE ACEL(AS ,AR,A W,DRAG,LIFT,THR, ALPHA, AOB.MASS) 
IMPLICIT REAL*8(A-HJ.-Z) 

C 

C COMPUTE THE ACCELERATIONS ON THE FLIGHT VEHICLE USING 
C EQUATIONS 3.28, 3.29, AND 3.30 
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c 

AS=(THR*COS(ALPHA)-DRAG)/MASS 

AW=(LIFT+THR*SIN(ALPHA))*SIN(AOB)/MASS 

AR=(LIFT+THR*SIN(ALPHA))*COS(AOB)/MASS 

C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C THIS SUBROUTINE IS THE COLLECTION OF DIFFERENTIAL EQUATIONS 
C WHICH DESCRIBE THE MOTION OF THE SPACECRAFT EQUATIONS 
C 3.9, 3.10, 3.1 1, 3.31,3.32, AND 3.33 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE ORB(X,XDOT,AS,AR,AW,G) 

IMPLICIT REAL*8(A-H0.-Z) 

DIMENSION X(6),XDOT(6) 

C 

XDOT(l )=X(4)*SIN(X(5)) 

XDOT(2)=X(4)*COS(X(5))*COS(X(6))/(X(l)*COS(X(3))) 

XDOT(3)=X(4)*COS(X(5))*SIN(X(6))/X(l) 

XDOT(4)=AS-G*SIN(X(5)) 

XDOT(5)=(AR-G*COS(X(5))+X(4)*X(4)*COS(X(5))/X(l))/X(4) 

DUM=TAN(X(3))/X(1) 

XDOT(6)=AW/(COS(X(5))*X(4))-X(4)*COS(X(5))*COS(X(6))*DUM 

C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C THIS SUBROUTINE TAKES CARE OF ALL HARDCOPY OUTPUT 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

SUBROUTINE WRT(X,T,THR,MASS,ALPHA,AOB,RHO,N,M) 

IMPLICIT REAL*8(A-HT^-Z) 

DIMENSION X(6),XDOT(6) 

PI=3.14159265359D+0 

RADIUS=X(1)/1000. 

VELOCITY=X(4)/1000. 

GAMMA=X(5)*180./PI 
ALFA=ALPHA* 1 80./PI 
HEAT=RHO**N*X(4)**M*9.652D-5 
SIGMA=AOB*180./PI 
OINC=ACOSD(COS(X(3))*COS(X(6))) 
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DELTA=0.0 

IF (OINC .NE. 0.) DELTA=ASIND(TAN(X(3))yTAND(OINC)) 

OMEGA=X(2)* 1 80./PI-DELTA 

WRriE(l 3, 1 )T,RADIUS,VELOCITY,MASS,GAMMA,OINC,THR,HEAT, ALFA 
.SIGMA 

1 F0RMAT(1X,F6.UX,F8.3,1X,F6.4,2X,F7.2,1X,F5.3,1X,F6.3,1X,F7.1, 

*1X,E9.3,1X,F7.3,1X,F7.3) 

C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

THIS SUBROUTINE ATTACH A HEADER TO THE OUTPUTS SO I CAN KEEP 
TRACK OF THE DIFFERENT CHANGES 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

SUBROUTINE HDR(R,V,MASS,THR, GAMMA, CSE) 

IMPLICIT REAL*8(A-HJ.-Z) 

CHARACTER* 10 CSE 
DO 1 1=13,14 
IF (LEQ.14)GO TO 1 
WRITE(I,*) 

WRITEa,*)’ SELECTED INITIAL INPUT DATA:’ 

WRITE(I,*) 

WRITE(I,2)CSE 

WRITE(I,3)R,V,MASS,THR,GAMMA 

WRITEd,*) 

WRITE(I,*) 

1 CONTINUE 

2 FORMAT(1X,A10) 

3 FORMATC INITIAL RADIUS (METERS) ’,F10.2y,’ INITIAL 
♦VELOCITY (METERS/SEC) ’J^.2y,’ INITIAL MASS (KG) 

* 'T^.2,/,' INITIAL THRUST (NEWTONS) 
initial FLIGHT PATH ANGLE (RADS) ’,F7.5) 

C 

WRITE(13,*)’ TIME RADIUS VELOCITY MASS GAMMA INCLI 

* THRUST QDOT ALPHA AOB' 

WRITE(13,*)' (SEC) (KM) (KM/SEC) (KG) (DEG) (DEG) 

* (N) (J/M2S) (DEG) (DEG)' 

C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C THIS SUBROUTINE IS A FOURTH ORDER RUNGE-KUTTA ROUTINE 
C TAKEN FROM PROF. ROSS SUBROUTINE 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE RK4(T,X,XDOT,NE,H,INDEX) 

IMPLICIT REAL*8(A-H0.-Z) 

INTEGER INDEX,! 

DIMENSION X(6),XDOT(6),SAVED(6),SAVEX(6) 

C 

INDEX=INDEX+1 
GO TO (1,2.3,4),INDEX 

1 DO10I=l,NE 
SAVEX(I)=Xa) 

SAVED(D=XDOT(I) 

10 X(I)=SAVEXa)+.5DO*H*XDOTa) 

T=T+.5D0*H 

RETURN 

C 

2 DO 20 1=1, NE 
SAVED(I)=SAVED(I)+2.DO*XDOT(I) 

20 X(D=SAVEX(I)+.5DO*H*XDOTa) 

RETURN 

C 

3 DO 30 1=1, NE 
SAVED(I)=SAVED(I)+2.DO*XDOT(I) 

30 Xa)=SAVEX(I)+H*XDOT(I) 

T=T+.5D0*H 

RETURN 

C 

4 DO40I=l,NE 

40 Xa)=SAVEX(I)+H/6.DO*(SAVEDa)+XDOT(I)) 

INDEX=0 

RETURN 

END 
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SAMPLE INPUT DATA FILE "DATA.DAT" 



This data file is for the orbital parameters of the chosen flight vehicle The name of the 
variable is on the left and the value placed on the right using a D13.6 format 



TYPE 

RADIUS METERS 
THETA RADIANS 
PHI RADIANS 
VELOCITY M/S 
FLIGHT PATH ANGLE 
PSI RADIANS 
BEGIN TIME SEC 
END TIME SEC 


AEROBANG 

6.4517440D+06 

O.OOOOOOOD+00 

O.OOOOOOOD+00 

7.6922000D+03 

O.OOOOOOOD+00 

O.OOOOOOOD+00 

O.OOOOOOOD+00 

6.0000000D+02 



TIME INTERVAL SEC 1 .OOOOOOOD-02 
PRINT TIME INTERVAL l.OOOOOOOD+00 
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SAMPLE INPUT FILE "AERO.DAT" 



This data file is for the aerodynamic parameters of the chosen flight 
vehicle The name of the variable is on the left and the value placed on 
the right using a D 13.6 format. Coefficients 1-6 are the lift and drag 
binomial coefficient 



COEFnCffiNTCl 
COEFFICffiNT C2 
COEFnCffiNT C3 
COEFnCffiNTC4 
COEFnCffiNT C5 
COEFnCffiNT C6 
AOA FOR CL/CD MAX 
AOB FOR BANG MANU 
DENSITY EXPONENT 
VELOCITY EXPONENT 
REFERENCE AREA 
INITIAL MASS (KG) 
FINAL MASS (KG) 
SPECEFIC IMPULSE (S) 
THRUST (N) 



1.500000D-02 

9.850000D-01 

9.240000D-01 

9.700000D-02 

-3.340000D-01 

2.619000D+00 

2.440000D-01 

1.570796D+00 

0.500000D+00 

3.150000D+00 

1.648000D+01 

5.185000D+03 

5.081000D+03 

2.950000D+02 

1.469700D+04 
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SAMPLE OUTPUT FILE "OUT.DAT” 



SELECTED INITIAL INPUT DATA: 



AEROBANG 

INITIAL RADIUS (METERS) 6445000.00 

INITIAL VELOCITY (METERS/SEC) 7621.30 

INITIAL MASS (KG) 5185.00 

INITIAL THRUST (NEWTONS) 14697.00 

INITIAL FLIGHT PATH ANGLE (RADS) 0.00170 



TIME RADIUS VELOCITY MASS GAMMA INCLI THRUST QDOT ALPHA AOB 



(SEC) 


(KM) 


(KM/SEC) (KG) 


(DEG) 


0.0 


6445.000 


7.6213 


5185.00 


0.097 


2.0 


6445.025 


7.6255 


5174.85 


0.089 


4.0 


6445.047 


7.6294 


5164.69 


0.080 


6.0 


6445.067 


7.6328 


5154.53 


0.072 


8.0 


6445.085 


7.6359 


5144.37 


0.063 


10.0 


6445.101 


7.6386 


5134.20 


0.055 


12.0 


6445.115 


7.6409 


5124.04 


0.047 


14.0 


6445.126 


7.6429 


5113.88 


0.039 


16.0 


6445.136 


7.6445 


5103.72 


0.031 


18.0 


6445.143 


7.6457 


5093.56 


0.023 


20.0 


6445.148 


7.6466 


5083.40 


0.015 


22.0 


6445.151 


7.6471 


5073.24 


0.008 


24.0 


6445.152 


7.6473 


5063.08 


0.000 


26.0 


6445.151 


7.6471 


5052.91 


-.008 


28.0 


6445.148 


7.6465 


5042.75 


-.016 


30.0 


6445.143 


7.6456 


5032.59 


-.024 


32.0 


6445.135 


7.6444 


5022.43 


-.032 


34.0 


6445.126 


7.6427 


5012.27 


-.040 



(DEG) (N) (J/M2S) (DEG) (DEG) 

0.000 14697.0 0.142E+07 3.602 90.000 
0.022 14697.0 0.142E+07 9.203 90.000 
0.053 14697.0 0.142E+07 11.580 90.000 
0.091 14697.0 0.142E+07 13.385 90.000 
0.135 14697.0 0.142E+07 14.892 90.000 
0.183 14697.0 0.142E+07 16.207 90.000 
0.236 14697.0 0.142E+07 17.386 90.000 
0.292 14697.0 0.142E+07 18.461 90.000 
0.353 14697.0 0.142E+07 19.454 90.000 
0.417 14697.0 0.142E+07 20.380 90.000 
0.484 14697.0 0.142E+07 21.251 90.000 
0.555 14697.0 0.142E+07 22.074 90.000 
0.629 14697.0 0.142E+07 22.857 90.000 
0.706 14697.0 0.142E+07 23.604 90.000 
0.786 14697.0 0.142E+07 24.321 90.000 
0.869 14697.0 0.142E+07 25.010 90.000 
0.955 14697.0 0.142E+07 25.675 90.000 
1.044 14697.0 0.142E+07 26.318 90.000 
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